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Optimum-Interval Punched-Card Tables 


This paper describes a new principle in table-making that leads to greater 
efficiency in punched-card applications. Tables of the ordinary kind, pro- 
ceeding by equal intervals of the argument, are often much more extensive 
than is necessary. In the new tables the number of entries is reduced to the 
minimum by allowing the interval of the argument to vary and by modify- 
ing the tabulated quantities so that interpolation is still performed in the 
usual way. 

The use of tables in computations that are made by means of punched 
cards usually entails the following operations:! the detail cards are sorted 
into the table cards according to the left-hand part of the argument, the 
function and first difference are gang-punched from the table cards onto the 
detail cards, and the interpolation of the function to the value corresponding 
to the complete argument is accomplished by one operation of the form 
Fand, with the multiplying punch. This assumes that the table is suffi- 
ciently extensive so that second differences are negligible. 

In some cases a less extensive table, used with an interpolation formula 
that includes second differences, may be more efficient. The necessity of 
deriving coefficients of the second difference may be avoided if the differ- 
ences in the table are replaced by the coefficients of the powers of m in the 
interpolation formula. Thus Bessel’s formula becomes 


n(n — 1) 
4 


F = Fy nh,’ 
= Fy + n(D, + nD.), 


where D; = A,’ — + Ay’), De = + A;”’) are the quantities to 
be tabulated instead of the customary first and second differences. This 
principle can easily be extended to higher orders of differences. 

The method just described is most efficient in application when the 
second differences are large and when the third differences may be neglected. 
When the contribution of the second differences becomes small it becomes 
relatively wasteful to take account of it in this way. Sometimes the diffi- 
culty may be overcome by some approximate method of sorting and gang- 
punching, but this is time-consuming and subject to error. The ordinary 
alternative is to extend the argument one more decimal place and to increase 
the number of cards in the table tenfold, a procedure which is even more 
wasteful and inefficient. 

Our purpose here is to describe a much more powerful expedient, the 
optimum-interval method. The interval of the argument is allowed to vary 
throughout the table, the criterion being that it shall be just small enough 
at every point so that linear interpolation will be legitimate, with the result 
that the number of cards in the table is greatly reduced. This advantage 
is obtained without complicating the use of the table in any way, and in 
spite of the unusual tabular intervals. 

In using any ordinary table the interpolating factor, m, is obtained by 
subtracting the tabular argument from the given argument, and by dividing 
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(Ao” + Ay’) 
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this difference by the tabular interval of the argument. When the tabular 
interval is an integral power of 10, then m is simply those figures of the 
complete argument beyond the position corresponding to the last digit of 
the tabulated argument. When the interval of the argument is not an inte- 
gral power of 10, let us suppose all the differences of the table to be divided 
by the interval of the argument. Form the product of this modified differ- 
ence by those figures of the complete argument beyond and including the 
position corresponding to the last digit of the tabulated argument. If this 
product is added to the function, the result differs from the correct inter- 
polate by an amount equal to the product of the modified difference times 
the last figure of the tabulated argument. By modifying each function in 
the table to the extent of deducting this product, the resulting table will 
admit of direct interpolation in the usual manner in spite of the odd interval; 
it is only necessary to retain one more figure in the multiplier. 

The value of the interval, w, is determined by the equation A” = w*d?f/dx?, 
where A” = 4, if the neglected second difference effect is not to exceed one- 
half unit. The sequence of intervals must be chosen so that every argument 
the last figure of which is zero is retained in the table. Otherwise there will 
be cases in which the full amount by which the function has been modified 
will not be restored. Of all the possible combinations of intervals, the follow- 
ing four (in addition to the unit interval) are all that it is necessary to 
employ in order to obtain the advantages of the method: 2, 2, 2, 2, 2; 
2, 3, 3, 3; 3, 3, 4; S, 5. 

A table of this kind differs from the ordinary table in respect of the 
accuracy attainable with a given number of figures. If the quantities in an 
ordinary table are accurate to 0.5 unit and if the second differences are 
negligible it is well known that a rounded-off interpolate has a maximum 
error of 1.0. The differences in the modified table are multiplied by factors 
larger than unity and the error of the interpolate depends on the interval of 
tabulation. The maximum rounding-off error is (1 + /g) units, where p 
is the largest interval of the argument and g is the unit of the argument. 
In many applications an error of this size is not important but if it can not 
be tolerated the remedy is to carry the values of the modified functions and 
differences to one more significant figure. In punched-card applications the 
additional figure usually does not complicate the use of the table or reduce 
the efficiency. 

To illustrate the advantages of optimum-interval tables we have pre- 
pared a six-place table of reciprocals, portions of which are shown below. 
The table is actually carried to seven places in order to insure six-place 
accuracy in the interpolates, and the tabular intervals have been chosen so 
as to keep the maximum error due to neglect of the second difference less 
than 5 units of the seventh place. Since the function is x1, A” = 2w*x-* = 4 
or x* = w’/2. Taking w = 2, 3, 4, 5, 10, 20, 30, 40, in succession, the corre- 
sponding values of the argument at which these intervals may be employed 
are, with sufficient approximation, 1.260, 1.660, 2.000, 2.320, 3.700, 6.000, 
7.700, 9.300, respectively. Let us consider a given argument lying between 
1264 and 1266. In the unmodified table it is necessary first to subtract 1264 
from the given argument and to divide by 2 in order to obtain the inter- 
polating factor n. In the modified table the differences are already divided 
by the interval. If the end figures of the argument beyond 126 are used as 
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an interpolating factor the result will be in error by 4D, and this amount 
has been subtracted algebraically from the true reciprocal of 1264 in order 
to obtain the function of the modified table. Thus the reciprocal of 12655 
is 7936393 — 5.5 X 6249 = 790202. It is advisable in practice to add 5 to 
every entry of the table in order to provide for automatic rounding-off of 
the interpolates to six figures and this has been done here. 

The total number of cards in the table of reciprocals is 1368, as compared 
with 9000 in the usual table, thus affording a saving of 2} hours of sorting 
and gang-punching each time the table is used. The loss due to an extra 
figure in the multiplier is one hour for each 10,000 detail cards. In the range 
from 1000 to 3700 and 370 to 600, the multiplier contains all the figures of 
the argument beyond the third figure, and from 600 to 999 it contains all 
beyond the second. In the range from 370 to 600, the modified table does 
not differ from an ordinary table. 

Another example of the improved efficiency possible with the use of 
optimum-interval tables is a table of seven-place sines and cosines. The 
usual table, with a tabular interval of 0°01, contains 9000 cards. The modi- 
fied table contains 2700 cards, with intervals 0°03, 0°03, 0°04. The maximum 
contribution of the neglected second difference is 0.6 unit of the seventh 
decimal. Each time this table is used, more than two hours of sorting and 
gang-punching time is saved. 

The application of this principle will be found of value in all cases of 
punched-card tables that differ considerably from linearity over the whole 
range of the argument. 


Arg. Unmodified Table Modified Table 

= 1/x f -D 
1250 8000005 6395 8000005 6395 
1251 7993610 6385 7999995 6385 
1252 7987225 6374 7999973 6374 
1253 7980851 6364 7999943 6364 
1254 7974487 6355 7999907 6355 
1255 7968132 6344 7999852 6344 
1256 7961788 6334 7999792 6334 
1257 7955454 6323 7999715 6323 
1258 7949131 6314 7999643 6314 
1259 7942817 6304 7999553 6304 
1260 7936513 12578 7936513 6289 
1262 7923935 12538 7936473 6269 
1264 7911397 12498 7936393 6249 
1266 7898899 12459 7936279 6230 
1268 7886440 12419 7936120 6210 
1270 7874021 12381 7874021 6191 
1990 5025131 5046 5025131 2523 
1992 5020085 7549 5025117 2516 
1995 5012536 5020 5025086 2510 
1997 5007516 7511 5025044 2504 
2000 5000005 7489 5000005 2496 
2003 4992516 7466 4999983 2489 
2006 4985050 9921 4999930 2480 
2010 4975129 7414 4975129 2471 
3690 2710032 3667 2710032 * 733 
3695 2706365 3657 2710020 731 


370 2702708 7285 2702708 7285 


' 
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Arg. Unmodified Table Modified Table 

x 1/x -A’ -D 
371 2695423 7246 2695423 7246 
597 1675047 2801 1675047 2801 
598 1672246 2792 1672246 2792 
599 1669454 2782 1669454 2782 
600 1666672 5537 1666672 2769 
602 1661135 5501 1666637 2751 
604 1655634 5464 1666562 2732 
606 1650170 5428 1666454 2714 
608 1644742 5393 1666318 2697 
610 1639349 5357 1639349 2679 


Pau HERGET and G. M. CLEMENCE 
U. S. Naval Observatory, 
Washington, D. C. 


1 Compare W. J. EcKERtT, Punched Card Methods in Scientific Computation, New York, 
The Thomas J. Watson Astronomical Computing Bureau, Columbia University, 1940. 
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138[A].— MATHEMATICAL TABLES Project, New York, A. N. Lown, tech- 
nical director, Table of Reciprocals of the Integers from 100,000 through 
200,009. New York, Columbia University Press, 1943, viii, 201 p. 
19.6 X 26.5 cm. Reproduced by a photo offset process. $4.00. 


This work is designed to provide a 7-place table of reciprocals between 100,000 and 
200,000, which will ‘‘expand by tenfold the scope of the existing tables in this interval.” 

The table is presented in what James Henderson has called the ‘‘modern arrangement”’ 
where the first five figures of the argument proceed in natural order down the page and the 
final figures from 0 to 9 are given along the top. Only the significant figures in the reciprocals 
are recorded. Since the entries are close together the differences vary from 100 at the 
beginning of the table to 25 at the end. Tables of proportional parts are given at the bottom 
of each page. 

According to the preface we learn that: ‘Preparation of manuscript tables was begun 
in December, 1934, by Dr. C. C. Kiess of the National Bureau of Standards. The work was 
completed under his direction in May, 1939, with the co-operation, at various times, of 
Messrs. B. F. Scribner, H. R. Mullin, W. G. Esmond, and J. Waldron.” 

It may be well to add a short account of the development of tables of reciprocals, which 
have been in the making for more than a century. The first adequate table appears to 
have been that of PETER BARLow, published in the first edition of his Tables in 1814 
[MTAC, p. 17], where the reciprocals are given for the first 10,000 integers. A. De Morgan, 
in the edition of 1840, says: ‘‘I cannot ascertain that any tables of square roots, cube roots 
or reciprocals comparable in extent to those of Mr. Barlow were ever printed before his.” 
In the fourth edition, issued in 1941 by L. J. Comrie, the reciprocals were extended from 
10,000 to 12,500. Comrie found ‘60 errors of a unit in the last decimal, but none greater” 
in the 1840 edition. 

The next significant contribution was R. PicarTE’s, La division réduite @ une Addition, 
Paris, 1861, which provided values of the reciprocals from 1000 to 10,000 to 10 significant 
figures, together with the first nine multiples of them. 

W. H. Oakes published his 7-place Table of the Reciprocals of Numbers from 1 to 100,000, 
in London in 1865. The volume under review is a continuation of this work, which was 
similarly printed in the ‘‘modern arrangement,” with differences at the side of the page. 
The more recent table of M. B. Corswortu, erroneously entitled ‘“‘Reciprocals for All 
Numbers from 1 to 10,000,000,” ! is also a 7-figure table over the same range as that of Oakes. 
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J. C. Houzeau in Acad. r. des Sciences . . . de Belgique, Bull., v. 40, 1875, p. 107, 
published a 20-place table of reciprocals of the first 100 integers and their first 9 multiples 
to 12 places. Apparently ignorant of this work, Kart PEarson, in v. 2 of his Tables for 
Statisticians and Biometricians, London, 1931, p. 257 gave 20-decimal values of the recip- 
rocals to 100. Ten-place reciprocals of the first 1,000 numbers are also given by E. G&LIn in 
his Recueil de tables numériques, Huy, 1895. 

Mention has already been made in this journal (p. 164) of the 17-figure tables for the 
reciprocals of the first 1,000 integers which are in the possession of the reviewer. He has 
recently acquired a 9-figure table of the reciprocals of numbers from 1 to 10,000, the work 
of CHarLEs MiLts of Northwestern University. 

1 A misstatement in the book under review. A more correct title of the work in question 
is as follows: Cotsworth’s Reciprocals for All Numbers From 1 to 10,000,000 with Direct Index 
Reference to Seven Significant Figures, complete for Numbers under 100,000 beyond which the 


extension to 10,000,000 zs rendered easy by the Additive Differences printed at the foot of each 
Table, York, England, M. B. Cotsworth, 1902.—Epirors. 


139[C].—H. S. Unter, “Natural logarithms of small prime numbers,” Nat. 
Acad. Sci., Proc., v. 29, 1943, p. 319-325. 17.1 XK 25.6 cm. 


The results here given are supplementary to those in Mr. Uhler’s publications already 
reviewed in MTAC, p. 55, 56, 20 (RMT 95, 96, 86). It was thus already noted that J. C. 
Apams (1887) calculated each of the quantities In2, In3, In5, and In7, to 272D; and that 
Mr. Uhler extended each of these to 328D and found In17 to 331D, and In71 to 213D. In 
the present paper the following new results are presented: In11 (329D), In13 (290D), 
In19 (299D), In23 (295D), 1n29 (300D), In31 (292D), In37 (291D), In101 (329D), and 
1n9901 (303D). Furthermore, each of the following is given to 155D and is supplementary 
to T. 4 of the monograph on 137-place values of In(i--n-10-*), RMT 86: in41, In43, 1n59, 
1n61, 1n67, In73, In79, 1n83, and 1n89. Thus for every N, prime and < 100, In N is given to 
at least 155D. The descriptions of checks applied tend to inspire confidence in the accuracy 
of the results. The author wished us to add, however, that in 1n43, line 3, 12917, should read 
12971; a slip due to failure of the printer to make this correction, indicated in proof; there 
is no error in the reprints. 


140(C, L].—WaLTER MEIssNER, Tafel der In T'-Function mit komplexem Ar- 
gumentbereich, Dresden Diss., Deutsche Mathematik, v. 4, 1939, p. 537- 
555. 20.8 X 28.9 cm. 

For the 192 complex arguments 
2 = (2m + i 2nv3)/24, and z = [(2m + 1) + 2(2m + 1)v3]/24, 
with m = 6,7, ---, 16,17; = 0, 1, ---, 6, 7 (triangular net), and for the gamma function 
T(z) = evtiv = 10-e%, 

u and v are given to 7D (except for 24 entries to 15D); \ = u loge is given to 7D, and 

¢ = v180/z is expressed in degrees, minutes, and whole seconds. (See F. EMpE, Z. angew. 

Math. Mech., v. 20, 1940, p. 295). The tables were intended as a supplement to the little 

tables of '(z) and In F(z) (p. 111, 109), z real, in the chapter on gamma functions in P. E. 


Boéumer, Differenzengleichungen und bestimmte Integrale, Leipzig, Kohler, 1939, where the 
complex case is discussed. A, 


141[D].—J. C. P. MiLLer, Tables for Converting Rectangular to Polar Co- 
ordinates. London, Scientific Computing Service, 23 Bedford Square, 
London W.C. 1, 1939, 16 p. 15.3 X 25 cm. Two shillings. Authorized 
American reprint, Dover Publications, 31 East 27th St., New York 16 
[1943]. Seventy-five cents. 
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The main object of these excellent tables is to facilitate the conversion of rectangular 
coordinates (x,y) to polar coordinates (r,@) by means of the relations r = (x? + 4°), 
= (y/x). 

Denote by / the larger and by s the smaller of |x| and |y|, and evaluate k = s/I; 
this is the argument of the tables and evidently lies in the range 0 to 1. Then we have 
r= (P+ =1(1+ and = tan*k =tan'(s/l), = cot 'k = cot (s/l) 
= 90° — 6°. Then @ = tan™ (y/x) may be determined in terms of 69 or ¢o if the result is 
required in degrees or in terms of 69 alone if radians are used, with the help of an “‘octant”’ 
scheme, in which the signs and relative size of x and y constitute the argument. 

The small interval of tabulation, namely 0.001, leads to the advantage that the dif- 
ference between two successive values of any function never exceeds 10; interpolation may 
therefore be performed mentally. For cases where r is required to three significant figures 
only, and @ to 0°.1 only, interpolation is unnecessary. 

The text discusses ‘‘Method with a calculating machine,”’ ‘Slide rule method” and 
“Applications of the tables.’’ Under the last heading it is noted that the tables may be used 
for transformation of harmonic constants a and b, obtained by harmonic analysis, to am- 
plitude c and phase angle e, in accordance with the relation 


csin (nt + €) = asin nt + bcos nt. 


The use of the tables for the evaluation of impedance and phase lag or lead from resistance 
and reactance will immediately suggest itself to electrical engineers. The evaluation of the 
magnitude and direction of a vector from rectangular components and the conversion of 
complex numbers from the form x + ty to the form r = e*@ are further applications. 


142[D, E]—MATHEMATICAL TABLES Project, New York, A. N. Lowan, 
technical director, Table of Circular and Hyperbolic Tangents and Co- 
tangents for Radian Arguments. New York, Columbia University Press, 
1943, xxxviii, 410 p. 19.6 X 26.5 cm. Reproduced by a photo offset 
process. $5.00. 


The main table of this interesting volume is a 400-page table of tan x, cot x, tanh x, 
coth x, for x = [0.0000(0.0001)2.000] radians. The number of decimal places varies from 5 
to 13 for the different functions and the different argument ranges. The greater part of the 
table gives values to 8 significant figures. Near poles more significant figures are given as 
noted below. There are auxiliary tables as follows: Table II gives these functions for 
x = [0.0(0.1)10.0 radians; 10D], conversion tables: degrees, minutes, and seconds to radians 
and vice versa,! integral multiples of +/2 to 15D. Interpolation coefficients p(1 — p)/2 
and p(1 — #*)/6 for p = [0.000(0.001) 1.000; 6D}. 

Except near x = 0 and x = 7/2, the main table was constructed, one entry at a time, 
by actual division of the corresponding values of sin x, cos x, sinh x, cosh x, taken from 
this Project’s table of circular and hyperbolic sines and cosines.? The values thus obtained 
were then verified by the reciprocal relations 


tan x cot x = tanh xcothx = 1. 


The resulting table was then differenced, to reveal errors in copying into the manuscript. 
This meticulous procedure, characteristic of this Project’s fine large tables, was made pos- 
sible by the large man- and machine-power available. Although the computer, to whom a 
good computing machine is available, can calculate isolated entries in the same way, or 
better still for the hyperbolic functions, by use of the formulas 


tanh x = (e* — 1)/(e* + 1), coth x = (e* + 1)/(e* — 1), 


where the value of e** can be taken from this Project’s excellent table of the exponential 
function,’ nevertheless the present table is more convenient (especiaily for the circular 
functions) when some interpolation is required. 
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The arrangement of the function values is a little confusing to the reader who merely 
thumbs the pages. The last three digits of each entry are separated from the earlier digits 
by a space. This makes for a variable number of digits between this space and the decimal, 
an irritating and perhaps hazardous circumstance. For x = [0.0000(0.0001)0.1000], all four 
functions are given to 8 decimals. Since cot x and coth x are quite large here, these functions 
are given with great accuracy, and since they are decreasing rapidly, interpolation is dif- 
ficult. For the sub-range x = [0.0300(0.0001)0.1000] a column of second differences is 
provided ranging from 74075 to 2000 units of the 8th decimal place. For still smaller values 
of x direct interpolation is impractical. The rest of the main table gives the four functions 
to 8 significant figures except for page 316 (x = [1.5700(0.0001)1.5750]) where cot x is 
given to 13 decimal places. On p. viii of the foreword, by H. T. D., one reads “‘The figures 
seen on the printed page may appear dull reading to the uninitiated. . . .” For once, this 
is not true of p. 316. Here one sees the function tan x passing through its agonizing crisis 
at x = 2/2, its onetime companion tanh x moving slowly ahead with apparent unconcern. 
Striking as these values of tan x are, they are of little practical value. To use any one of 
them, the value of x must be known with extreme certainty. Interpolation is not possible. 
Perhaps it would have been more practical to tabulate, near x/2, the function 


tan x — (x/2 — x)" 
which is interpolable. Similarly near x = 0 the functions 
cot x — 1/x, coth x — 1/x, 


might have been tabulated. The policy of basing the table on 8 significant figures rather 
than a fixed number of decimals, is presumably to make uniform the relative error, rather 
than the actual error in the tabular entries. This expedient fails to produce the desired 
result in the case of tanh x and coth x. Near the end of the table the tanh x values are, 
on the average, nearly ten times as accurate as those of coth x, on the basis either of relative 
or actual error (since both functions are nearly unity). For the last half of the table, the 
relative error in tanh x averages about sixteen percent that of coth x. 

The introduction contains (p. xxiii-xxxviii) a bibliography of tables and charts of 
circular and hyperbolic tangents and cotangents and their logarithms, the inverse hyper- 
bolic functions, the gudermannian and related functions, together with tables of roots of 
equations involving tan x. No short bibliography could cover the very large number of 
titles of this description. In the present case only 73 important items are listed. These do 
not include any tables produced by New York Tables Project. Of these 73 tables only 17 
refer to circular or hyperbolic tangents or cotangents given to more than 6 decimal places. 
Of these 17, only one table is for radian measure. This is the unreliable table of Hayasut* 
in which tan x and tanh x are tabulated up tox = 50, the main table being for x = [0.0010- 
(0.0001)0.100(0.001)3.00(0.01)9.99] Thus the table under review fills a real gap in the biblio- 
graphy of these elementary functions. It will, no doubt, be useful chiefly in problems in 
which a large number of values of one of the four functions are needed quickly. The fact 
that linear interpolation gives such good results over most of the table makes it quite useful 
to the engineer or physicist. It will help to wean away the typical astronomer and civil 
engineer from his sexagesimal table of log tan x. 


D.H.L. 


1In the reviewer's opinion such conversion tables should be used only for pencil and 
paper work. With a machine available they are too slow and dangerous to bother with. 

2 Tables of Circular and Hyperbolic Sines and Cosines, New York, 1939. 

3 Tables of the Exponential Function e*, New York, 1939. 

*K. Hayasut, Sieben- und ——— Tafeln der Kreis- und Hyperbelfunktionen und 
deren Produkte sowie der Gammafunktion, Berlin, Springer, 1926. 


Nore.—R. C. A. requests that his statement concerning works of Ives and BENSON in re- 
views of MTAC, 1943, p. 8-10, be considered rather than the one quoted on p. xxx, I. 16-20, 
from Scripta Mathematica, v. 2, 1933, p. 91. 
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143[D, E, L].—HERBErT Bristot Dwicut (1885-— ), Mathematical Tables 
of Elementary and Some Higher Mathematical Functions including Trigo- 
nometric Functions of Decimals of Degrees and Logarithms, New York 
and London, McGraw-Hill Book Co., 1941. Third impression (with addi- 
tions), 1944, viii, 231 + 8 p. 15.2 X 22.8 cm. $2.50. Reproduced by a 
photo offset process. 


A detailed statement of the contents of this volume will best exhibit its value. First 
differences for interpolation are given in connection with practically all of the tables. 
After the tables are many references to sources where more extended tables may be con- 
sulted. Tables I-V: sine, cosine, tangent, cotangent and their logarithms, and also secant 
and cosecant for every hundredth of a degree, to 5D, p. 1-101. T. VI-VIII: sine, cosine, 
for each 0°.001, 0 to 2, and tangent, 0 to 1.570, p. 102-113; T. IX—XI: sin™ x, cos x, 
tan“! x, x = [0°.000(0°.001) 17.000; 4D] and x for more extended values, p. 114-121; 
T. XII: Inx, x = [1.000(0.001)3.00(0.01)10; 4D], p. 122-127; T. XIII-XIV, e*, e-, 
x = [0.000(0.001)5; 46S], p. 128-147; T. XV-XVII: sinh x, cosh x, tanh x, x = [0.000- 
(0.001)3; 4-6S], p. 148-165; T. XVIII-XIX: x, cosh! x, x = [0.000(0.001)2.00(0.01)- 
10; 4-5S], p. 166-175; T. XX: tanh x, x = [0.000(0.001)1; 4—5S], p. 176-178; T. X XI: 
binomial coefficients 1C, to 25C,, p. 179; T. XXII: (a? + 6*)4/a, b/a = 0.000(0.001)1, with 
a reference to Miller’s Tables (see RMT 141), p. 180-181; T. XXIII: factorials !, n = 1 
to 50, 1-6S, p. 182; T. XXIV: log m!, m = 1 to 250, 5-7D, p. 182-183; T. XXV: Gregory- 
Newton interpolation coefficients, »C,, r = 2(1)6, p = [0.01(0.01)1.00; 2-7D], p. 184-185; 
Lagrangean interpolation coefficients for interpolating without columns of differences, 
p. 186-187; T. XXVI-XXVII: surface zonal harmonics P,(x), Pn(cos 6), m = 1(1)10, 
x = [0.00(0.01)1; 1-6D], @ = 0°(1°)90°; p. 188-195; T. XXVIII: first derivatives of the 
surface zonal harmonics P,(cos 0), @ = 0°(1°)90° from Farr’s article, and therefore very 
erroneous (see MTE 33), p. 196-197; T. XXIX: complete elliptic integral of the first kind, 
p. 199-203; T. XXX: complete elliptic integral of the second kind, p. 204-205; T. XXXI: 
Bernoulli’s numbers (35), p. 206; T. XXXII: Euler’s numbers (35), p. 207; T. XX XIII: 
I'(n), m = [1.00(0.01)2; 5D], p. 208-209; T. XXXIV: erf x, x = [0.000(0.001)1.50(0.01)3; 
5D], p. 212-213; T. XX XV: ber x, bei x, ber’ x, bei’ x, p. 214-217; T. XXXVI: ber x + 
ibei x, ber’ x + ibei’ x, p. 218-221; T. XXXVII-XXXIX: Riemann zeta function ¢(s), 
s = [— 24.0(0-1)24; 10-11S], also tables of (s — 1)¢(s) with &, 5¢ and ¢’(s)/¢(s), p. 222- 
227; T. XX XX: log x, x = [1.00(0.01)9.99; 4D], p. 228-229; Index, p. 231. The zeta func- 
tion tables were furnished by H.T.D. 

Such are the contents of the first edition. In the current edition Preface and Contents 
are printed from type rather than “‘offset’’; literature notes have been modified on pages 
195, 197 and 213; on p. 124 In 2-269, for -8193 has been substituted .8194; argument 
headings on p. 208-209 have been corrected; and 8 pages have been added. On the new 
p. 195A-B are the first derivatives of surface zonal harmonics P’,(x), m = 2(1)8, 
x = (0.00(0.01)1.00; 2-6S]; the erroneous p. 196-197 persist. Six new pages, 213A-F, are 
devoted to a table of the normal probability integral. 


144[F]—Hansray Gupta, “On the class-numbers of binary quadratic 
forms,’’ Tucuman, Argentina, Universidad, Revista, s. A, Matem. y Fisica 
Teorica, v. 3, 1942, p. 283-299. 17.8 X 27.1 cm. 


An integral binary quadratic form is a homogeneous expression 
(1) (x1, X2) = ax + bxix2 + = [a, c] 


with integra! coefficients a, b,c. If the greatest common divisor (a, b,c) of a, b, and c is 
equal to unity, then f is said to be properly primitive. A binary quadratic form F(y1, y2) 
obtained from f(x1, x2) by a substitution x; = airy: + aizye (i = 1, 2) with integral coef- 
ficients and of determinant a11022 — a@2ia12 = 1, is said to be properly equivalent to f. Both 
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f and F have the same discriminant b* — 4ac = — d. All properly equivalent forms are said 
to form a class. 

The present paper contains a table of the number h’(d) of the classes of properly 
primitive binary quadratic forms of negative discriminant — d, 1 < d < 12500, together 
with the description of the procedure employed in the computation. It is the most extensive 
class-number table known to the reviewer. However, it contains no such related information 
on genera as do, for example, the extensive tables of Gauss.! The total number h(d) of 
classes of discriminant — d is obtained from the well-known formula 


h(d) = = h'(d/F). 


There exist closed expressions? of theoretical importance for h’(d). None of these, how- 
ever, are really suitable for computation. The method of construction of the present table 
is based on the following considerations. As is well known, every properly primitive quad- 
ratic form of negative discriminant — d is equivalent to one and only one reduced form f 
for which*® 


(2) -a<b€a O<a&c, b>O0ifa=c, (b,c =1. 


Thus the number h’(d) is the same as the number of properly primitive reduced forms of 
discriminant — d. Should one know the number‘ D(d, B) of reduced forms (2) with |b| = B, 
one may use the relation 
Wd) =ZDd,B), BS 
to obtain h’(d). 
To compute D(d, B) it suffices to determine the number N(i, B) of decompositions 
4 = (d + b*)/4 = a-c into two positive factors a and c such that 


(3)  (C(a,c),b) =1, 
where C(a, c) = (a,c). 

For, D(d, 0) = N(d/4, 0). Next let B # 0. 

In case i is not a perfect square, then 


(4) D(d, B) = 2N(i, B) + 1, if B|i and C(B, c) = (B,i/B) = 1, 
and 
(5) D(d, B) = 2N(i, B) otherwise. 


In case 7 is a perfect square, then formulae (4) and (5) hold if (B, i) > 1. If, however, 
(B, 1) = 1, the values in (4) and (5) must be increased by + 1. Thus, in this case: 


(6) D(d, 1)=2N(, B) + 2; Did, B) = 2N +1 if B >1 and either (B,i) = 1 or 
Bi and (B,i/B) = 1; D(d, B) = 2N for other values of B. 


These formulae follow at once from (2). The term + 1 in (4), accounts for the reduced 
form with a = B = b which is properly primitive only if C(a,c) = 1. The unit to be added 
when is a perfect square accounts for the form with a = c, b > 0. This form is properly 
primitive if (B, a) = (B, a?) = (B, 7) = 1. Actual computation may be arranged systemati- 
cally as follows. First comes the “preliminary table.” 


C(a, c) 1 2 4 1 4 2 
211282625 002020 2 O80 

Cc 1 3 5 1 15 

225 N 22.8 4.6.6 4.83 & & 


Here, in view of (3), N(i, B) is obtained by counting the number of C(a,c) with 
B <a < i, which are prime to B. The value of D is obtained through the formulae (4)-(6). 
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Next, the D entries are copied out opposite the proper values of d in another table and 
are totalled by rows to give h’(d). 

Thus, the method of tabulation used by the author is essentially the usual method of 
construction of the table of reduced forms in which the scheme of computation is so arranged 
that only the presence of a reduced form [a, b, c] is recorded and not the form itself. One 
observes with interest that the desirable scheme in this case differs from the excellent 
arrangement by Wright® in tabulating reduced forms themselves. As in Wright the method 
calls for only a simple succession of elementary operations. 

St. Louis University ARNOLD E. Ross 


1C. F. Gauss, Werke, v. 2, p. 450-476. For reviews of this and other similar tables 
nag pty LEHMER, Guide to Tables in the Theory of Numbers, Nat. Res. Council, Bull. 
no. 

2 Compare L. E. Dickson, History of the Theory of Numbers, v. 3, Washington, D. C., 
1923 (New York reprint, 1934), chapter VI. 

3 Compare L. E. Dickson, Modern Elementary Theory of Numbers, Chicago, 1939. 

4 The author uses D((d+ B*) /4, B). 

5 The author omits this row and records the values of D directly. 

®H. N. Wricut, “On a tabulation of reduced binary quadratic forms of a negative 
as,” California, University, Publs. in Math., v. 1, no. 5, 1914, p. 98-114+2 folding 
plates. 


145[F].—Cuao Ko and S. C. Wana, “Table of primitive positive quater- 
nary quadratic forms with determinants < 25,’’ Academia Sinica, Sci- 
ence Record, v. 1, nos. 1-2, Aug., 1942, p. 54-58, Chungking, China. 
[In the 260 pages of this publication 97 pages are devoted to 25 papers 
announcing results of mathematical research.] 18 X 25.4 cm. 


n 
In the study of diophantine equations (1) f(x) = m, where f(x) = ¥ aj,xix; is a 


quadratic form in m variables x1, ---, x, with integral coefficients a;;, it was noticed at an 

early date that integral solutions of (1) may be obtained from the solutions of a related 

equation (2) F(y) = m, where F(y) is a quadratic form equivalent to f(x), i.e., a form obtained 
n 


from f by a transformation 7: x; = > 7;,y; with integral coefficients 7;; and of determinant 
im1 


|r«j]| equal to unity. Thus, seeking solutions of (1) and (2) independently would involve 
needless duplication of effort which would be avoided if one would select for study only one 
of the unlimited number of equations (1) whose left-hand members f are equivalent. A 
process of selection of a limited number of representative forms, desirably a single one, for 
each set (class) of equivalent forms is called reduction, and the selected representatives are 
referred to as reduced forms. Since the number of classes of quadratic forms of a fixed 
determinant |a;;| is finite, a limited number of reduced forms suffices to represent all 
forms of a given determinant. 

In the paper under review, the authors have constructed a table of primitive (g c d of 
2a,; and a;; is unity) reduced positive quaternary quadratic forms of determinant D < 25, 
using a method of reduction which gave exactly one representative for each class of equiva- 
lent forms. Such a table should be considered an improvement of CHARVE’s! original table, 
since his method of reduction led to more than one reduced form for some classes. The 
authors were apparently unaware, however, of the existence of a similar table by S. B. 
Townes?, D < 25. Townes’s method of reduction employs ideas of Eisenstein and also leads 
to a unique representative. Reduced forms in the two tables differ only in the choice of the 
cross-product coefficients. 

There are a number of errata in thearticle under review; on p. 56f, D = 20, for (1,1,2,20), 
read (1,1,2,10); D = 21, for (12330,—10), read (12340, —10); and D = 24, for 
(2283 —10; —10), read (2233 —10; —10). On p. 58, the fourth determinant is 9 
and not 8. In addition to these errata, noted by the authors in an errata list at the end of 
the issue of the Record, we note that the entries (2234 —1:0; —1) for D = 16 and 


de 
Pp. 


conm 


(2 
of 
clz 
th 
in 
+ 
sh 
| i 
x 
n 
4 I 
| 
( 
| 
t 
I 
i 
= 
> ! 

| 


RECENT MATHEMATICAL TABLES 183 


(22330 —120 —12) and (2235 —1:0; —1) for D = 24 are incorrect, being actually 
of determinants 28, 11 and 36 respectively. A representative for one of the twenty-four 
classes of determinant 24 is missing, and the entry (2 33 3 —1 0, —1 02) of D = 25, should 
read (2233 —1 02). 

To restore completeness of the table one may use the following reduced forms from 
Townes in place of the entries referred to above: the form (2 23 3 1.013) for D = 16, and 
the forms (2 23412013), (2333 —12101:), and (12440; 2) for D = 24. 

The count of the number of classes yields 14 for D = 16 and 24 for D = 24. We have, 
in the above, used the well-known notation (abcdrstuvw) = ax;? + bx? + cx;? + dx? 
+ 2rxexs + 2sxixs + + Zuxixg + + A run of zeros is indicated by 
0, for brevity. One may note that the entry (1 1440; 2) for D = 16 in Townes’s table 
should read (1 1450, 2). 

St. Louis University ARNOLD E. Ross 


1L. Cuarve, ‘Table des formes quadratiques quaternaires positives réduites dont le 
joe est égal ou inférieur 4 20,’ Acad. d. Sci., Paris, Comptes Rendus, v. 96, 1883, 
p. 773-775. 

2S. B. Townes, ‘‘Table of reduced positive quaternary quadratic forms,”’ Annals Math., 
s. 2, v. 41, 1940, p. 57-58. 
EpritoriAt Note: Mr. Ross was formerly called ARNOLD CHAIMOVICH, as we learn from his 
mimeographed English translation from the Russian of P.S. Nazrmov’s Prize Essay (Mos- 
cow, ae ~ pplications of the Theory of Elliptic Functions to the Theory of Numbers, Evans- 
ton, Ill., 1928. 


146[F].—D. H. Leumer, ‘‘Ramanuja.’s function r(m),”" Duke Math. J., v. 
10, 1943, p. 483-492. 


The numerical function 7r(m) is the coefficient of x* in the power series expansion of 
x[(1 — x)(1 — x*)(1 — x4) ---].24 This function has long been of interest in the theory of 
numbers, particularly in the arithmetical applications of the elliptic modular functions. 
If m, m are coprime, 


(1) 1(mn) =1(m)r(n); 
if p is prime 
(2) r(p**!) = r(p)r(p*) — pUr(p™), a=1. 


Hence r(n) is calculable from the values of 7(p*), which in turn are calculable for a = 2, 3, 

- when 7(p) is known. An explicit formula for 7(p) is lacking. Ramanujan conjectured 
(from scanty numerical evidence, namely, the 10 cases determined by p < 30) that 
|r(p)| < 2p"/2, so that = 2 cos 6p, 6, real. This unproved conjecture, whose 
truth was regarded by Ramanujan as “highly probable,” is called by Hardy the “Ramanujan 
hypothesis.”’ The author’s purpose is best described in his own words: “In seeking to dis- 
prove the Ramanujan hypothesis I have examined all primes p < 300, that is, the first 
46 primes, as well as p = 571, and I find that in all these cases the hypothesis is true. There 
is only one ‘near miss’ which occurs at p = 103.” 

The table gives the value of r(m) for m = 300. For composite m the values were com- 
puted by (1), (2). For prime a recurrence formula in which the arguments of r(m) decrease 
by successive pentagonal numbers is used. The formula is obtained by an application of 
Euler’s power series expansion of II(1 — x"), m = 1, 2, --- to the logarithmic derivative of 
the generating identity for r(), and is 

(nm — 1)r(m) = D[m — 1 — §.25(3m? + m)| [n — 4(3m* + m)], 
the summation referring to 1 = |m| < an, 6a, = 1 + (1 + 24n)!. This formula has prac- 
tical computational advantages over Ramanujan’s recurrence formula involving triangular 
numbers, derived similarly from Jacobi’s expansion of II(1 — x*)’. 
The calculations for r(m), m prime, were checked by Ramanujan’s congruence 


mod 691, 
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where o1:(”) is the sum of the 11th powers of the divisors of . The table as a whole was 
checked by comparison at frequent intervals of the values modulo 691 of Er(s), Zou(s), 
s =x. Further congruences, to be discussed in another paper, also were used as checks. 
In addition to the table and the method of computing and checking it, the paper discusses 
problems connected with the orders of r(m), Zr(s), neither of which is known. 


E. T. BELL 
California Institute of Technology 


147[I]—ARTHUR SHEPARD LITTLE (1872- ), A Table of Interpolation 
Multipliers, for obtaining through the means of calculating machines, Inter- 
mediate Rates Bond Values at yield intervals of one ten-thousandth percent. 
Boston, Financial Publishing Co., [1927], 29 p. 17.5 X 21.6 cm. $4.00. 


The main table of this work gives the four coefficients A, B, C, D, in the following 
special case of Lagrange’s interpolation formula: 


ai S(@ + 2hx) = A(x)f(a) + B(x)f(a + h) — C(x)f(a + 2h) + D(x)f(a + 3h) 
w 

A(x) =— (2x — 1)(% — 1)(2x — 3)/3, 

B(x) = 2x(x — 1)(2x — 3), 

C(x) = x(2x — 1)(2x — 3), 

D(x) = 2x(2x — 1)(x — 1)/3, 


for x = [0.000(0.001)0.500; 9D]. The reason for the peculiar ‘‘2x’’ in the above is that the 
table is intended for use with bond tables having an argument interval of 0.05 percent. 
There is also on p. 29 a small table giving the three coefficients in the corresponding three- 
term formula: 
f(a + 2hx) = A(x)f(a) + B(x)f(a + h) — C(x)f(a + 2h) 

where 

A(x) = (2x — 1)(@ — 1), 

B(x) =— 4x(x — 1), 

C(x) =— x(2x — 1), 


for x = [0.00(0.01)0.50; 4D]. The author’s 8 pages of introduction are interesting (not to 
say amusing) but inadequate. None of the above formulas appears in the book and not 
even one small example is worked out to show the business man how to use these tables. 
Instead there are repeated warnings against failing to notice that formulas (which are not 
given) contain minus signs so that the crank of the calculating machine must be turned in 
the negative direction at the proper time. To make this point very clear the coefficients 
C(x) in the two tables are printed in red and the column headings bear the injunction 
‘*Turn Off,’’ while the other columns are headed ‘Turn On.” 

The author does not indicate how he discovered these interpolation formulas. There is 
no reference to Lagrange or to any other source. The author believes that his method is 
“absolutely original and is now being made public for the first time.’’ 

The Lagrange type of interpolation formula is, I believe, destined to become more and 
more used by computers to whom a good calculating machine is available. The more tra- 
ditional formulas employing differences are products of an age of pencil and paper cal- 
culation. Many tables are printed without differences. When the successive differences are 
available or worked out, the various terms of the interpolation formula vary so in size 
that they are a constant source of irritation if not of actual error. With the Lagrange type, 
the separate terms are of the same order of magnitude and are automatically accumulated 
in the product register. A small table of Lagrangean interpolation coefficients, such as 
the one under review,! should be in the library of every computer. 

D. H. L. 

1 Perhaps one of the best elementary introductions to the subject may be found in 
K. Pearson’s Tracts for Computers, nos. II-III, Cambridge, University Press, 1920, with 
illustration examples and short tables of coefficients in the 4 to 11-point ‘“‘Lagrangean 
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formulae.” LAGRANGE gave the formula now associated with his name in Paris, Ecole 
Polytechnique, J., v. 2, 1795, p. 274-278. But this formula was discovered twice earlier, 
namely: (1) by Epwarp Warinc, R. So. London, Trans., v. 69, 1779, p. 59-67, and also 
in generalized forms; (2) by L. EULER, Opuscula Analytica, v. 1, 1783, p. 184, setting forth 
the general Lagrangean formula for equal intervals. 

Reference may be given also to T. L. Kettey, The Kelley Statistical Tables, New 
York, Macmillan, 1938. There are cubic (4-point, [0.000(0.001)1; 10D]), quintic (6-point, 
[0.00(0.01)1; 10D]), and septic eae, [0.0(0.1)9; 11D]) interpolation tables. See RMT 
130. And finally we may note MATHEMATICAL TABLES Project, Table of Lagrangean 
Coefficients, Ordnance Department, 1941, 40 p. 2226.5 cm. This edition of a 5-point 
table [0.000(0.001)2; 7D] was not available to the public until prepared for distribution 
7 the Marchant Co., Oakland, California, in 1942. This is now published on p. 383-403 
of MATHEMATICAL TABLES Project, Table of Bessel Functions Jo(z) and J;(z) for Complex 
Arguments, New York, Columbia University Press, 1943; see RMT 151. The Project has 
a the = _ elaborate volume, entitled Tables of Lagrangian Interpolation Coefficients 

XXxii, p.). 


Eprror1aL Note.—We include here a review of this table published before 1933, for two 
reasons, because only the slightest attention has been paid to it in any mathematical 
periodical, and because there is little doubt that the author, working most of his life in a 
clerical capacity for the First National Bank and one of its affiliates, in St. Louis, was an 
independent discoverer of the utility of these so-called Lagrange in tion coefficients. 
He published at least three other volumes (1915-26), listed in the Library of Congress 
Catalogue, and a good many articles in Journal of Accountancy. 


148[I].—A. J. THompson, Tables of the Coefficients of Everett's Central- 
Difference Interpolation Formula. Tracts for Computers, No. V. Edited 
by E. S. Pearson and published by the Department of Statistics, Uni- 
versity College, London, Cambridge University Press, second ed., 1943. 
viii, 32 p. 16 X 23.3 cm. 5 shillings. See also MTE 41. 


The first edition of Mr. Thompson’s tables was published in 1921 and was quickly sold 
out. The present edition differs considerably from the first (xvi, 21 p.). The size of the 
introduction is cut by half and its contents (somewhat paradoxically) are increased. Also 
there are five tables of coefficients in the new edition, the first three columns of Table I 
covering the whole of the single table in the original edition. The reputation of the dis- 
tinguished author guarantees the accuracy of the figures published. It may be expected, 
therefore, that the tract will be in great demand. 

A considerable portion of the original Introduction was devoted to numerical examples 
illustrating in varying degrees the same type of application of Everett’s interpolation for- 
mula and coefficients. In the second edition the number of numerical examples is reduced 
to just a single one. In addition to this, the author gives to the reader the benefit of his 
wide experience in dealing with various difficult situations. A straight application of 
Everett’s formula is possible when the table, in which it is desired to interpolate, contains 
the central differences. This, of course, is far from being always the case. A convenient 
method is indicated by which explicit use of the central differences can be avoided. 

Other material in the new Introduction deals with the construction of mathematical 
tables. The first step suggested is to compute directly the values of the function to be 
tabulated for a basic framework of widely spaced values of the argument. These computa- 
tions should yield some three or four places of decimals in excess of the accuracy of the 
proposed table. Once the framework is completed, the use of Tables II, III, [V and V of the 
new edition of Thompson’s tract will permit a relatively easy process for filling of the gaps 
by reducing the original increment of the argument in the ratio of one to ten or more. The 
coefficients of Everett’s formula are used to compute both the values of the function and 
the corresponding central differences. 

Table I gives the values of the first four coefficients for values of the argument 0 < @ < 1 
proceeding by increments of .001. It also gives the central differences of the coefficients 
which could be used for interpolation among the values published. Table II contains coef- 
ficients of the same orders as Table I but proceeding at intervals of @ equal to .01. Table 
III gives the values of six coefficients, ¢2 to €12, with complete differences, but proceeds at 
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increments of @ equal to .1. Table IV has the same order of coefficients and the same incre- 
ment of @ as Table III but the range of 0 is from —5 to +6. The primary purpose of this 
table is to facilitate interpolation in other tables near the ends of the range of their argu- 
ments, where the central differences of higher orders cannot be made available without 
previously extending the range of the original table. In Table V the coefficients extend to 
€is and proceed at increments of @ equal to .2. 4 

Those who have frequent occasions to interpolate in tables in which linear or quadratic 
interpolation is inadequate, and have had to do so without the benefit of Mr. Thompson’s 
tract, may think that its title is not satisfactory and that it should be changed to “Inter- 
polation Made a Pleasure,” or to ‘‘A Hard Life Made Easy.” 


J. NEYMAN 
University of California, Berkeley 


149[I, J] —H. E. Sauzer, “Coefficients for numerical differentiation with 
central differences,” J. Math. Phys., M.1.T., v. 22, 1943, p. 115-135. 
17.5 X 25.5 cm. 


This paper provides the coefficients for the expansion of the mth derivative of ¢(x), 
when ¢(x) is expanded in terms of central differences, that is to say, for the coefficients of 
and in the expansions: 

vok 


von 


= + Re, 


where R; and Rz are remainders, whose explicit forms are found in standard treatises on 
finite differences. For the mean differences of odd order, we replace the customary notation, 
yé”-! or 116-1, by the simpler notation 6”—!. The original expansion of ¢(x) in terms of 
central differences is usually referred to as the Newton-Stirling series. 

The author defines the coefficients in terms of B,™(x), the Bernoulli polynomial of 
mth order and mth degree, as follows: 


Ge — 2n)! and Ay = — — Qn)! Bg” ..,(»). 


It is unfortunate that he did not also state them in terms of the familiar differential 
coefficients of zero, namely, 


= — 1)!, and = Dol") /(2y)!. 


2n—1 


Coefficients of even order can be immediately written down from those of odd order 
by the obvious relationship 
Ax, = (m/v) . 


The author used the following new identity in a partial check of his computations: 


1 

In order to achieve his purpose in making “‘this table the ‘ultimate’ in coefficients for 
numerical differentiation,” the author computes the coefficients as far as the 52nd deriva- 
tive. For the first 30 derivatives the author gives the exact values, in common fractions, up 
to the difference of 30th order, and also for some coefficients of differences beyond the 30th. 
“For all derivatives beyond the 30th, exact values are given for coefficients of differences 
going as far as some difference between the 41st and 52nd. Elsewhere, that is, for most of 
the coefficients of the 31st to 42nd differences, 18 significant figures are given, with accuracy 
to within 0.6 unit in the last significant figure.” 
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To put it otherwise, the derivative formulas can now be written to at least 18 significant 
figures for the first 32 derivatives up to differences of 42nd order, the 33rd, 34th, 35th, 
36th, 37th, 38th, 39th, 40th, and 41st derivatives up to orders 43, 44, 45, 46, 47, 48, 49, 50, 
and 51 respectively, and the derivatives of orders from 42 to 52 to differences of 52nd order. 

The coefficients were checked to about 10 significant figures by means of the recurrence 
formula given above. 

Although it is difficult to see where these tables would be used in the calculation of 
derivatives of such high order, since differences beyond eight or ten are almost never 
encountered, nevertheless the coefficients themselves are related to some interesting func- 
tions and could quite conceivably be of importance in problems other than those for which 
the table has been computed. We are greatly indebted to the author for providing us with 
these new constants. 


150[L]—Car_ AuGust HEuMAN (1870- ), “Tables of complete elliptic 
integrals,” J. Math. Physics, M.1.T., v. 20, 1941, p. 127-206. 17.5 x 
25.4 cm. 


The first two tables in this paper give values to 6D of the complete elliptic integrals 
of the first and second kinds, Fo(a) and Eo(a), for a = 0°.0(0°.1)90°, with first differences; 
except that the differences are omitted for a > 65° in the case of Fo(a). The functions 
Fo(a) and Eo(a) are 2/x times the Legendre integrals F and EZ. Linear interpolation is 
sufficiently accurate except in the range a > 65° for Fo(a). To fill this gap there is given 
a table of an auxiliary function in which linear interpolation can be used: 


Gola) = Fo(a) + (2/z) In (90° — a®) 


for a = 65°.0(0°.1)90°, with first differences. 
An extensive table of complete elliptic integrals of the third kind Ao(a, 8) is also given, 
p. 160-197, where 
Ao(a, 8) = Fo(a)E(a’, B) — { Fo(a) — Eo(a)} F(a’, 8). 


Here a’ = x/2 — a and F(a, 8) and E(a, 8) are the incomplete integrals of the first and 
second kinds. Tabular values are for a = 0°(1°)90°, 8 = 0°(1°)90°, and for a = 0°.0(0°.1)5°.9, 
B = 80°(1°)89°, all to 6D. 

Heuman obtained his values of Fo(«) and Eo(a) directly from values of log F and 
log E in Table I of Legendre’s Exercices de calcul intégral sur divers ordres de transcendantes 
et sur les quadratures, Paris, 1816, Tome 3, and Traité des fonctions elliptiques, Paris, 1826, 
Tome 2, with the aid of Vega’s ten-place Thesaurus Logarithmorum. The values of Ao(a, 8) 
were computed from Legendre’s Table I, to 12D-14D, for each tenth of a degree, and 
Table IX, the latter being a table of the incomplete integrals F(a, 8) and E(a, 8) to 9D- 
10D, for each degree. A list of 42 errors in Legendre’s tables is given, most of them being 
in the incomplete integrals. 

An appendix to the paper contains applications to the motion of a spherical (or more 
general) pendulum, and the gyroscopic pendulum. In the former case a chart is given for 
the apsidal angle as function of the two independent variables 2; and z2, which are the 
greatest and least heights attained by the bob during its motion. The gyroscopic pendulum 
is considered in the case where the spindle starts at an inclination 4 with zero initial velocity. 
Two charts are given, for the apsidal and zonal angles respectively, as functions of the two 
independent variables 4 and yw where yu is the coefficient of stability. 


P. W. Ketcuum 


151[L]— MATHEMATICAL TABLES Project, New York, A. N. Lowan, tech- 
nical director, Table of the Bessel Functions Jo(z) and J;(z) for Complex 
Arguments, New York, Columbia University Press, 1943, xliv, 403 p. 
19.6 X 26.5 cm. Reproduced by a photo offset process. $5.00. 
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The main part of this useful volume consists of a table of real and imaginary parts 
of Jo(z) and J:(z) as functions of z = pe”, for p = [0.00(0.01)10; 10D], # = 0°(5°)90°. 
The last twenty pages contain a table of the five-point Lagrangean Interpolation Coef- 
ficients for p = [0.000(0.001)1; 10D], twice previously published (MTAC, p. 94). There 
is also a short foreword by H. BATEMAN, and a long Introduction by Mr. Lowan, followed 
by a Bibliography (65 items) of Tables of Bessel functions and of applications of Bessel 
functions of a complex argument of the form rit. The four pages of “‘Contour Lines” for 
Jo(z) = Jo(pe*) = U(p, 6) + iV(p, ¢), Jo(2) = Jo(pe*) = Re, Ji(z) = Ji(pe*) = u(p, ¢) 
+ iv(p, ¢), Ji(z) = Ji(pe*) = Se*’, are of considerable interest. The publication of the 
Bessel function tabulation makes possible the numerical solution of a large number of 
problems, which have hitherto only been discussed symbolically. Apart from tables with 
pure imaginary arguments, and for arguments in the form ri#, DinnrK’s short tables! 
of Jo(z) and J:(z) are the only ones previously published. 

Bessel functions are ubiquitous in mathematicai physics, from acoustics to stochastics. 
When the wave equation in more than one dimension is separated in rectangular coor- 
dinates, the separated equations have an irregular singular point at infinity, and the 
solutions are exponential or trigonometric functions. If the coordinate system involves a 
point center or a single axis of revolution, such as for spherical or circular cylindrical coor- 
dinates for instance, the radial equation has a regular singular point at the origin and an 
irregular singular point at infinity. The general solution of this type equation is the con- 
fluent hypergeometric function, but in the great majority of cases of practical interest, the 
special case of the Bessel equation is the result. There are other more complicated coor- 
dinate systems in which the wave equation separates, but the resulting ordinary differential 
equations have a more complicated system of singular points, and the solutions (Mathieu 
functions, Lamé functions and the like) are not well enough tabulated to make them 
readily usable in specific problems. 

The peculiar position of the Bessel function is again evident when one studies the 
applications of the Laplace transform, whether in wave or diffusion problems or in prob- 
ability theory. In many coordinate systems the characteristic integral expression 


fe cos 


occurs. This can be shown to be proportional to one or another of the solutions of the 
Bessel equation of order m, depending on the choice of the contour for the integral. 

Perhaps part of the reason for the preeminence of the Bessel function in mathematical 
physics is that it is neither too complex nor too simple. A solution of the wave equation 
expressed in terms of trigonometric functions would not be adequate to use in studying 
the behavior of waves from point or line sources, which are best expressed in terms of 
Bessel functions. On the other hand, although a solution of the oblate spheroidal waves 
from a circular disk would be even more valuable than a solution for a point source, the 
oblate spheroidal functions needed for such a solution are so much more difficult to com- 
pute than are the Bessel functions, that no adequate tables are yet available from which 
to obtain complete numerical results. Bessel function solutions are complex enough to 
exhibit certain general wave properties, and yet simple enough to tempt the computer. 

As the Bibliography well suggests, tables of Bessel functions for real and for pure 
imaginary arguments are now fairly numerous. A great number of problems can be solved 
numerically by their use, and the tables under review do not represent an advance in this 
direction. There are a number of important problems, however, whose solution requires 
values of the Bessel functions for complex values of the argument; and in this field the 
present Table represents a distinct advance. 

Uses for the functions, for the argument evi, are fairly well known (the ber, bei and 
related functions), and a number of tables have been published of these functions.? Problems 
relating to the behavior of electromagnetic waves in cylindrical conductors (skin-effect) 
require the use of these functions, as also do calculations of the fluctuation of temperature 
in a cylinder whose surface temperature is a sinusoidal function of time (furnace with 
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thermostatic control). The present tables are more complete than any of the earlier ones 
of ber x, bei x (the B.A.A.S. Tables, for instance) for this special case of # = 45°. 

For other values of 3, however, the present tables are unique, and their publication 
makes possible calculations of considerable immediate interest. The behavior of micro- 
waves in cylinders of dielectric material in which there is energy dissipation is one appli- 
cation; and the scattering of electromagnetic waves from such a cylinder is another. The 
attenuation of acoustic waves down hollow cylinders* can now be calculated, and computing 
the transmission of sound in a shallow sea with absorbing bottom is made considerably 
easier. In fact, wherever cylindrical waves cause energy dissipation, either in the medium 
or at the boundary, Bessel functions of complex argument are encountered, and the present 
tables are of use. The companion volume, giving values of No and N; for complex values 
of the argument, will also be most useful when it is published. 

The reviewer has not had the leisure to test the accuracy of the present Tables. How- 
ever, the accuracy of the previous publications of the Mathematical Tables Project is well 
known. 


Pattie M. Morse 
Massachusetts Institute of Technology 


1A. N. Dinntk, All Russian Central Committee, [Communications], 1922, p. 121-126; 
in K. Havasu, Fiinfstellige Funktionentafeln ... , Berlin, Springer, 1930, p. 
105-109. 

2 A.G. WEBsTER, B.A.A.S., Report, 1912, p. 56-68; ber x and bei x and their derivatives, 
x =[0.0(0.1)10.0; 9D], with 7 differences. For corrections by H. G. SavinGeE see B.A.A.S., 
Report, 1916, p. 122; there are also corrections by H. B. Dwicur. SavipGE gave tables of 
ker, kei and their first derivatives and other tables in B.A.A.S., Reports, 1915, p. 36-38, 
and 1916, p. 108-121. An abridgement of some of the tables in these Reports is given by 
H. B. Dwicat in (a) Tables of Integrals and Other Mathematical Data, New York, Mac- 
millan, 1934; see RMT 154; and (b) Mathematical Tables, New York, McGraw-Hill, 1941; 
see RMT 143. The latter has 5-place tables of functions and derivatives for # = 45°, 135°. 

Morse, ‘The transmission of sound inside pipes,”” Acoustical So. Am., J., v. 11, 
1939, p. 205. 


152[L].—E. O. Powe t, “An integral related to the radiation integrals,” 
Phil. Mag., s. 7, v. 34, Sept. 1943, p. 600-607. 17 X 25.3 cm. 


The following function is tabulated 


Inydy/(y 1) = (— — 1)? 2), 


= InxIn (1 — x) — 2x" (0 < x 1); 
1 


Ri(x) + RI(1/x) = 1/2(log x)?. From these relations the values of Ri(x) were computed 
for x = [0.00(0.01)2.00(0.02)6.00; 7D], with second central differences, negative throughout. 
The author writes as follows: ‘It is hoped that the tabular values will be found correct to 
six decimal places; the seventh should not usually be in error by more than one unit.’’ The 
“differences are sufficient to give seven-figure accuracy when x is greater than 0.2, six- 
figure between 0.1 and 0.2; below this limit they increase very rapidly.” 

The integral /i7 In ydy/(y — 1) was recently encountered in a physical problem.* 
This integral is equal to Réi_,[In (1/x)] — 2°/6, where Rin(x) = /7°dy/y"(ev — 1), an 
integral tabulated by Arrey for 6 values of m (see MTAC, p. 140, no. 46). DeByE gave?* 
a short table of (1/x)So*ydy/(e¥ — 1) = 2*/(6x) — (1/x)Ri_s(x). 

1S. R. Finn and E. O. Pow t, ‘‘The chemical and physical investigation of germicidal 
aerosols. II: The aerosol centrifuge,’ J. Hygiene, v. 42, 1942, p. 364. 


2P,. Desye, “Interferenz von Réntgenstrahlen und Warmebewegung,” Annalen d. 
Phys., s. 4, v. 43, 1914, p. 85-86. 
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153[L].—F. VANpDREY, ‘‘Tafel der acht ersten Kugelfunktionen zweiter Art,” 
Z. angew. Math. Mech., v. 20, 1940, p. 277-279. 20.8 X 29.5 cm. 


These tables emanated from the aerodynamic research laboratory at Géttingen. They 
present the values of the Legendre spherical functions of the second kind Q,(x), from 
Qo(x) to Qr(x), for x = [0.00(0.01)1.00; 5D]. By virtue of the relation Q,(— x) 
= (— 1)"*1Q,(x) the table covers also the range — 1 =x =0. 

In calculating the table the values of the function 


Qo(x) = tanh x = § In [(1 + x)/(1 — x)] 


were taken from the 9-place table of K. Havasu, Sieben- und mehrstellige Tafeln der Kreis- 
und Hyperbelfunktionen, Berlin, Springer, 1926, p. 9-87. Q:(x) to Q,(x) were calculated from 
the following equations: Qi(x) = xQo(x) — 1, Q(x) = P2(x)Qo(x) — (3/2)x, Qs(x) 
= Ps(x)Qo(x) — (S/2)x* + %, Qu(x) = Pa(x)Qo(x) — (35/8)x* + (55/24)x. Qs to Qr were 
then determined by means of the recurrence formula, 


mQn(x) + (m — 1)Qn-a(x) — (2m — 1)xQns(x) = 0. 


In the tables a dot placed after the upper part of the fifth figure indicates that this figure 
has been increased. Linear interpolation to 3 places is possible in Qo(x) to Qs(x) from 0 to 
about 0.86, and in Q;(x) from 0 to 0.69. 


154[L, M].—H. B. Dwicut, Tables of Integrals and other Mathematical Data, 
New York, Macmillan, 1934, x, 222 p. 14 X 21.5 cm. $1.75. See also 
RMT 30. 


The first important American table of integrals was the one by the brilliant Harvard 
physicist, B. O. Perce (1854-1914) which started in a small way as a 32-page pamphlet 
(1889), and was bound in with the second edition of W. E. ByERLy’s Elements of the Integral 
Calculus, Boston, 1889. After tremendous labor this, A Short Table of Integrals, was ex- 
panded to a book of 144 pages (Boston, 1910). The third edition (156 p.), revised by W. F. 
Oscoop, appeared in 1929. After 479 indefinite integrals, there are 44 miscellaneous definite 
integrals, followed by more than 60 formulae under the heading Elliptic Integrals. Then 
follow (p. 73-115, nos. 570-938) various auxiliary formulae in trigonometric, hyperbolic, 
elliptic, and Bessel functions; series; derivatives; Green’s theorem and allied formulae; 
table of mathematical constants; general formulae of integration; note on interpolation. 
The final pages (116-154) contain various tables including those of the probability integral, 
of elliptic integrals, of hyperbolic functions, and of I'(m). 

And now we are considering a volume by a professor of electrical engineering at Mas- 
sachusetts Institute of Technology. In a general way the volume of Dwight is not greatly 
dissimilar to that of Peirce but it is considerably more elaborate. Dwight has much more 
detailed numbering of formulae! (with many cross references), and careful statements of 
ranges for which they are valid. On p. 219-220 are listed 35 volumes of tables, treatises, 
etc., to which references are often given as authorities for statements made, or as sources 
where more elaborate results appear. Peirce’s tables in this list are naturally referred to 
several times. Tables 1050, of ber, bei, ber’, bei’, ker, kei, ker’, kei’, are taken from those of 
Webster and Savidge in B.A.A.S. Reports, 1912, p. 56-68; 1915, p. 36-38; 1916, p. 122. 
Dwight has corrected the value for bei 8.9 given by Webster as — 8.002--- instead of 
— 28.002---. Dwight made this correction first in Amer. Inst. Electr. Engs., Jn. v. 42, 
1923, p. 830. But he failed to correct Webster’s value for bei’ 3.7 0.134 686 760 which 
should have been 0.131 486 760 (Dwight, A.I.E.E., Trans., v. 58, 1939, p. 787). On pages 
215-217 of Table 1050 are original tables of ber, x, bein x, ber’n x, bei’, x, ker, x, kein x, 
ker’, x, kei’, x, for m = 1(1)5; bern x + i bei, x = Jn(xi!) = i"J,(xit); ker, x + i kein x 
=1-"K,,(xi}). The first four of these tables, 3D-6D, were first published in A.I.E.E., Trans., 
v. 42, 1923, p. 858, and reprinted in Trans., v. 48, 1929, p. 814-815; of the second four, 
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mostly 6D-8D, the greater part of the values were calculated by W. H. Hastincs as thesis 
work at M.I.T. in 1927. 

The section on Bessel Functions, nos. 800-845, is a very useful collection of recurrence 
formulae, identities (functions and integrals), with real and complex arguments. The 
notation of Gray and Mathews and Macrobert’s A Treatise on Bessel Functions, London, 
1922 (1931), is here employed. 

The section of formulae on Elliptic Functions and Integrals occupies nos. 750—789.2 
and in nos. 1040-1041 there are two brief tables of K and E. In K, 87°.6, for 4.562, read 4.561. 

The series and formulae nos. 1-50 include some dealing with Bernoulli and Euler 
numbers. 

This volume, by the author of several mathematical tables, is an excellent and exceed- 
ingly useful one, compiled with great care. The “Contents” and “Index” are wholly 
adequate. Some other errata are noted in MTE 32. 

A. 


1 We have already drawn attention (MTAC, p. 66) to a valuable volume, Smithsonian 
Mathematical Formulae and Tables of Elliptic Functions, 1922 (corrected reprint 1939), 
314 p. The Mathematical Formulae (p. 1-219) were prepared by the physicist Epwin P. 
Apams. He writes, “In order to keep the volume within reasonable bounds, no tables of 
indefinite and definite integrals have been included. For a brief collection, that of the late 
Professor B. O. Peirce can hardly be improved upon; and the elaborate collection of Bierens 
de Haan show how inadequate any brief tables of definite integrals would be. A short list 
of useful tables of this kind, as well as of other volumes, having an object similar to this 
one, is appended.” Nevertheless, there is much in common with the volumes under — 
Dwight doubtless received more than one suggestion from this volume which has a 
number of literature references. In the useful chapter on infinite series, p. 109-144, — is 
still one error, as Mr. W. D. Lambert, of the Coast and Geodetic Survey, has ‘recently 
pointed out; on p. 122, under 6.42, no. 4, the third and fifth terms of the right-hand member 
should each be preceded by the sign —. Legendre and Bessel functions are considered at 
some length on p. 191-219. 


155[L, M].—A. H. Heattey, “A short table of the Toronto function,” 
R. So. Canada, Trans., v. 37, sect. III, 1943, p. 13-29. 16.2 K 25 cm. 


By “Toronto function” the author designates the three-parameter function T(m, n, r) 
defined in terms of the four-parameter function 


T(m, n, p,a) = (2at)dt 
by means of the relation 
T(m, n, 7) = n, p, a), 


where r = a/p. Many relations between the Toronto function and other functions such as 
the confluent hypergeometric function M(a, y,x), Bessel functions of fractional orders, 
the error function, etc., are summarized in Table 1-A. The author gives also (without 
derivation) a number of recurrence formulae between three T’s whose parameters m and 
n differ by integers, as well as a number of formulae for special values of these parameters. 
He also gives the differential equation satisfied by T when considered as a function of r 
and an asymptotic expansion based on the confluent hypergeometric function. 

The short table included with the article gives 5-place values of T(m,n,1) for the 
following values of the parameters: 


m =— 3; =— 2(0.5)2; r = 0(0.2)2; 5, 6, 10, 25, 50 
m=0; n =— 2(0.5)2;r = 0(0.2)4; 5, 6, 10, 25, 50, 
m=%4; mn =— 2(0.5)2;r = 0(0.2)2; 4, 5, 6, 10, 25, 50, 


m=1; m =— 2(0.5)2;r = 0(0.1)2(0.2)3. 


The computations of T(m, n, r) were carried out on the basis of the relations in Table 
1-A for sufficiently small values of r and on the basis of the asymptotic expansion for suf- 
ficiently large r (mostly for r > 4). 


' 
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The computed values were checked with the aid of the recurrence formulae and oc- 
casional differencing. In order to compute T(m, n, r) for small values of r, the author found 
it necessary to prepare two small tables of the function M(a, y, x) and e*M(a, y, x). The 
first lists M(4, 3, x), M(4, 2, x), M(3,1,x), M(4,2,x) and M(%, 3,x) for 11 values of x 
ranging from 0 to 4 at irregular intervals. The second lists e~*M(j, 3, x), e-*M(, 2, x), 
e-*M(i, 1, x), e-?M(, 2, x) and e~*M( 3?, 3, x) for x ranging between 0 and 4 at the same 
intervals. There is a gap in the values of 7(— 4, 2,r) between r = 2 and r = 5 anda 
gap in the values of 7(4, 2, r) between r = 2 and r = 4, 

The reviewer wonders whether it would not have been more expedient to base the com- 
putation of T(m, n, r) for small values of r directly on the differential equation rather than 
on its expression in terms of the confluent hypergeometric functions whose series expansion 
is so slow that the author had to take 19 terms of the series. 

Conceivably the method developed by C. Lanczos (‘‘Trigonometric interpolation of 
empirical and analytical functions,” J. Math. Physics, M.I.T., v. 17, 1938, p. 123-199) 
might have yielded expressions valid in a range of r including the gaps above noted. 

ARNOLD N. Lowan 


156[L, S].—L. D. RosENBERG, “Metod rascheta zvukovykh polel, obrazo- 
vannykh raspredelennymi sistemami izluchatelei, rabotaiushchimi v 


zakrytykh pomeshcheniiakh”’ [Method of calculating sound fields gen- 
erated by distributed systems of radiators operating in closed places] 
Zhurnal Tekhnicheskoi Fiziki, v. 12, 1942, p. 247-248. 16.6 X 25.9 cm. 


These pages at the end of the article contain a table of T(u, a) = J,'**+%e-™™dx, 
for a = [0.0(0.1)0.7, 0.785; 5D], « = 0.0(0.1)1.0, 1.5, 2.0, 3.0, 3.5, 4, 5, 6, 8, 10, 15, 20, 30, 
40, 50, 70, 100. 


157[M].—A. N. Lowan and J. LADERMAN, ‘‘Table of Fourier coefficients,” 
J. Math. Phys., M.1.T., v. 22, 1943, pp. 136-147. 17.5 X 25.5 cm. 


In this paper we find the tabulated values of the two functions 
S(k,n) = f x* sin nxx dx, and C(k, n) = f x* cos nxx dx, 


over the range m = 1(1)100, and & = 0(1)10. The tables are given to 10 decimal places 
with a plus sign indicating that the eleventh decimal place is 5 or larger. This is a rather 
interesting innovation, since the usual rule is to provide in the last place a figure which is 
not in error by more than .5 units, the sign of the error being either plus or minus. In this 
connection it might be worth while calling attention to the device employed by L. M. 
MILNE-THOMSON and L. J. Comrie in their Standard Four-Figure Mathematical Tables, 
London, 1931, who say with respect to their tables that ‘‘a further indication of the figures 
omitted is given by placing a ‘high’ dot after the last figure if these omitted figures, in 
units of the last figure retained, lie between .1666... (i.e. $) and .5000..., and a ‘low’ dot 
if they lie between .5000... and .8333... (ie. 3).” 
The two functions tabulated satisfy the recurrence formulas: 


S(k, n) = S(1, 2) 


and the “cross relations”’: 


2m), = 48 62, ») coe 2, m); 


S(k, n) = SU, n) C(k — 1, n), C(k, n) =— = se — 1, #). 


As to the accuracy of the table we may quote the authors, who say: “The computations 
were performed originally to twenty decimal places with such extreme care, that the 
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‘cross relation’ checks on the twenty-place values did not reveal a single error. The present 
table is an abridged version of the original twenty-place manuscript.” 

The table is designed to facilitate the computation of Fourier coefficients for functions 
which may be represented by polynomials of degrees not exceeding 10. 

Since the table under review is unique in the field of harmonic analysis, which has been 
developed vigorously during the past forty years by astronomers, physicists, meteorologists, 
and economists, it may be of interest to indicate its relationship to other tables. 

In harmonic analysis the interest has focused upon the two sums, 


2nzxt 
N 


2nxt 


> fi cos and Ba Sesion 


or the corresponding integrals, 
2 2 
ft) cosas, and f $00) sin ay, 
a 0 a a ° a 


which are usually to be evaluated by finite integration. 
Hence most tables designed for use in harmonic analysis provide values of the functions 


M cos — and M sin —. 


The following is a bibliography of such tables: 


H. H. Turner, Tables for facilitating the use of harmonic analysis. London, Milford, 
1913, 46 p. 

L. Z1pPERER, Tafeln zur harmonischen Analyse periodischer Kurven, Berlin, Springer, 
1922, 12 p. 

L. W. PoLiak, Rechentafeln zur harmonischen Analyse, Leipzig, Barth, 1926, 23 + 
17 + 240 p. 

L. W. PoLiak, Handweiser zur harmonischen Analyse, Prague, Bursik & Kohout, 1928, 
72 + 98 p. 

P. TEREBEs!, Rechenschablonen fiir harmonische Analyse und Synthese. Berlin, Springer, 
1930, 13 + 11 p. + 12 tables. 

A. HussMann, Rechnerische Verfahren zur harmonischen Analyse und Synthese, mit 
Schablonen fiir eine Rechnung mit 12, 24, 36, oder 72 Ordinaten. Berlin, Springer, 1938. 28 p. 
+ 12 folding tables. 

K. Srumprr, Tafeln und Aufgaben zur harmonischen Analyse und Periodogramrechnung. 
Berlin, Springer, 1939, vii + 174 p. 

Since the demands of practical statistics do not usually require many decimal places, 
the above tables are computed mainly to three or four significant figures. Thus the last 
table, which is typical, contains three-place values of the two functions for M = 1, 
b = 3(2)39, nm = O(1)4p; p = 2(4)40, mn = O(1)ip; p = 4(4)40, m = O(1)ip; M = 0(1)1000, 
n = 1, p = 8, 12, 16, and 24; M = 2(2)100, p = 21(1)27, 29, 30(1)35, 37(1)39, m = 0(1)39; 
M = 1(1)100, p = 20(4)40, nm = 0(1)}p. 

Since occasionally greater accuracy is required than is provided by these tables the 
reviewer has prepared in manuscript form an 8-place table of the two functions for M = 1, 
b = 5(1)75, nm = O(1)p. 

It will be clear that the table which is the subject of this review will be useful in har- 
monic analysis in those instances when the statistical function f(#) can be approximated 
by a polynomial of degree not exceeding ten. The values of A’ and B’ are then immediately 
approximated by appropriate sums of the the tabulated values of C(k, m) and S(k, m), when 
we use the transformation t = ax. The table was computed during 1940-41 in connection 
with work of the MATHEMATICAL TABLES PROJECT. 


H. T. D. 
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References have been made to errata in RMT 143 (Dwicut), 145 (Ko 
and Wanc), 150 (LEGENDRE), 154 (ApAms, DwicHTt, WEBSTER), N18 
(SILBERSTEIN), N19 (CARRINGTON, JAHNKE & EmDE), and N20 (Boys, 
ByRNE, KNAPPEN). 


29. BESSEL and HANSEN Tables of Bessel Functions. 


It is well known that Bessel published in 1826 a table of J,(x), n = 0, 1; x = [0.00- 
(0.01)3.20; 10D]. Because of misinformation published in sources usually highly authori- 
tative it is not so well known that one of P. A. Hansen’s tables of 1843 was of Jn(x), n = 0, 
1; x = [0.0(0.1)20; 6D]. Among others this was reprinted by Lommel in 1868. 

In the absurd paragraph of 16 lines on ‘Tables de fonctions de Bessel” in Encyclopédie 
des Sciences Mathématiques, 1, 5, 2, p. 228, one reads (omitting footnote references) “F. W. 
Bessel a calculé des tables des fonctions Jo et J: pour des arguments variant de x = 0a 
x = 3,2. P. A. Hansen les a étendues jusqu’a x = 10, E. Lommel jusqu’a x = 20.”’ [Hansen 
did not extend Bessel’s 10-place table for hundredths of a unit to x = 10, and Lommel 
did not extend Hansen’s table.] In Watson’s Bessel Functions, 1922, p. 655, there are three 
misleading statements, ‘‘Hansen constructed a Table of Jo(x) and J:(x) to six places of 
decimals with a range from x = 0 to x = 10.0 with interval 0.1 . . . Hansen’s table was 
reprinted . . . by Lommel who extended it to x = 20.” E. M. Horsburgh repeats one 
misleading phrase and makes one new error, on p. 57 of his Modern Instruments and Methods 
of Calculation, 1914, ‘“‘P. A. Hansen’s extension of Bessel’s table is reproduced by . . . 
E. Lommel. . . . It gives Jo(x) and Ji(x) from x = 0 tox = 20 at intervals of .01 through- 
out the lower part of the range.” It would seem as if all three authors did not recognize 
that J,(x) = Ij, (Hansen) = J"(x) (Lommel) [= J? (Bessel)]. 

A. 


30. E. W. Brown with the assistance of H. B. HEprick, Tables of the Motion 
of the Moon, New Haven and London, 1919. Compare MTAC, p. 29. 


Apart from the erratum on p. xiv of these Tables, two lists of errata were published 
by Brown, namely: (1) in Astr. J., v. 34, 1922, p. 54; and (2) in Yale Univ., Astr. Obser- 
vatory, Trans., v. 3, 1926, p. 157; in 1. 4 from the bottom, for 92-l’, read 82-1’. Eight other 
errors sent to us are now presented. The one in Section II was supplied by D. H. SADLER, 
Superintendent of H. M. Nautical Almanac Office, three of those in Sections I and VI 
came from W. J. Eckert, Director of the Nautical Almanac Office, Washington, and the 
rest, in Sections I, III, VI, were taken by DirK Brouwer, Director of the Yale University 
Observatory, from two copies of the Tables used by the late E. W. Brown. 


Section I 
P. 44,1. 3, for 70, read 71. 
80, 1.-2, for n-, read n+, 
09, 1. 


1 10, for log cos w, read log sin w. 


Section II 
P. 31, arg. 84 for 1966, for 2305, read 2405. 
Section III 
P. 160, col. 17, arg. 16.5, for 19608, read 19598. 
P. 187, col. 153, arg. 9.0, for 23878, read 23868. 
Section VI 
P. 63, col. 47, arg. 180, 190, 200, for 110, 117, 112, read 210, 217, 212. 
P. 92, 0°59’0”", for 007 4464, read 007 4454. 


194 


Cc 
al 
fe 
ol 
fc 
h 
c 
n 
I 
: 
| 
\ 
} 


MATHEMATICAL TABLES—ERRATA 195 


31. Epwin CHAPPELL (1883-1938), A Table of Coefficients to Facilitate Inter- 
polation by Means of the Formulz of Gauss, Bessel and Everett. London, 
The Author, 1929, viii, 27 p. 22.3 K 28 cm. 


As the author of these tables is no longer living, it may be of interest to mention that 
the balance of the 110 copies printed (by Chappell himself) are now available from Scientific 
Computing Service Ltd., 23 Bedford Square, W.C. 1. 

Besides the two corrections given by the author on page 2, the formula on the second 
page of the Introduction for G** = B* should be preceded by a minus sign. On pages 26 
and 27 the values of G** = B* for 0.38 and 0.62 (which are the same) are both in error: 
for 0.004 56167, read 0.004 56157. The occurrence of this error (which was found about ten 
years ago by Mr. D. H. SapDLeEr) in two places shows that it is an error in recording the 
original calculations, rather than one of composition. The fact that the values in Table I 
for 0.380 and 0.620 are correct shows that a comparison between the two tables was not 
made. 

This table was compiled at the suggestion of the present writer, whose copy is in- 
scribed “3/110. With the author’s compliments and in recognition of the suggestion which 
has resulted in this book.”’ The calculations were done on the first Brunsviga Dupla cal- 
culating machine imported into England—a machine which, by curious fate, also came into 
my hands some years after Chappell’s death. 


Epitor1at Note: This well printed and bound little volume of Chappell is worthy of a 
place in any large mathematical library, along with his Five-figure Mathematical Tables: 
consisting of logs and cologs of numbers from 1 to 40,000, illogs (antilogs) of numbers from 
.0000 to .9999, lologs (logs of logs) of numbers from 0.00100 to 1,000, illologs (antilologs) of 
numbers from 6.0 to 0.5000, Together with an Explanatory Introduction and Numerous Ex- 
amples; also, trigonometrical functions and their logs of angles from 0°-90° at intervals of 
1 minute, With Subsidiary Tables. Edinburgh, London, Chambers, and New York, Van 
Nostrand, 1915, xvi, 320 p. Compare Q 4. Chappell was much interested in Pepysiana, and 
published (1933-36) at least four items concerning Samuel Pepys and his bibliography. 


32. H. B. Dwicut, Tables of Integrals and Other Mathematical Data, New 
York, 1934. See RMT 154. Most of the following errata were contributed 
by the author. 


P. vii, for page numbers 113, 114, 114, 117, 118, 121, read respectively: 111, 112, 112, 

115, 116, 119. 

93.1. Change 3aX to 3a?X. 

After the sub-headings on pages 101, 102 and 144, insert (@ > 0). 

Omit No. 591. 

603.2 For the part in square brackets, read [sin x # 0]; and for logsin x, read log|sin x|. 

603.4. For the part in square brackets, read [cos x ¥ 0]; and for logcos x, read log|cos x|. 
(See P. Franklin, Differential Equations, New York, 1933, p. 275). 

P. 127, last line, for x/4, read x. 

630.2, last line, read [x? < x*/4]. 

673.19, for 657.6 read 657.8. 

679.19, for 657.5, read 657.7. 

Since a statement is made in 810.8 about the size of the error from using the asymptotic 
series, it would be better to add (1/r)Ko(x) to 810.6; and (K,(x)/r)e***®* to 810.7. 

815.2 for s = 1, reads = 0. 

836.4 for ker x, read ker’x. 

850.1 for It, read The integral. 


- 851.4. Change to the asymptotic series: 


1 1 139 
+ 1) = E + 288n* 51840n* 


Omit No. 861.9. 
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Table 1030: read logio cosh .32 = .02187; sinh 1.20 = 1.5095; tanh 3.15 = .99633; tanh 3.55 
= .99835. 

Table 1050: read bei’ 3.7 = + 0.131 486 760. From x = 4.6 to 6.4, delete the last 3 digits 
of the values of bei’x. 

P. 214, last line, for 109, read 122. 


33. C. C. Farr, “On some expressions for the radial and axial components 
of the magnetic force in the interior of solenoids of circular cross-section,”’ 
R. So. London, Proc., v. 64, 1899, table p. 199-202. Reprinted without 
change in all five editions of JAHNKE and Empe, 1909-1943. 


On comparing with our recently calculated table for dP,(cos 6)/d0, m = 1(1)7, we found 
53 last figure errors, of which the 12 of more than a unit in the last figure are as follows: 
n = 3, 31°, for — 2.0654, read — 2.0656; n = 4, 5°, for — 0.8570, read — 0.8567; n = 4, 
64°, for 1.6306, read 1.6300; n = 5, 54°, for 2.0173, read 2.0178; n = 5, 64°, for 1.5420, 
read 1.5418; m = 5, 67°, for 1.1187, read 1.1183; n = 5, 83°, for — 1.4831, read — 1.4827; 
n = 6, 23°, for — 3.0917, read — 3.0902; n = 6, 48°, for 2.3715, read 2.3713; n = 6, 57°, 
for 1.1936, read 1.1934; m = 6, 83°, for — 1.4487, read — 1.4484; n = 7, 70°, for — 1.984, 
read — 1.934. None of these 53 errors occur in the corresponding table of H. Tallqvist, 
Grunderna af Teorin for Sferiska Funktioner jimte Anvandningar inom Fysiken, Helsingsfors, 
1905, p. 422-431. 
We checked also the table of P,(cos 6), m = 1(1)7, @ = 0°(1°)90°, on p. 121-123 of the 
fifth edition of JAHNKE and Empk, and found no error. 
MATHEMATICAL TABLES PROJECT 
New York City 


EprtroriaL Note.—The table of P,,(cos @) here referred to, is the one by HOLLAND, JONES 
and Lams, originally with 92 errors (MTAC, p. 136-137). These were faithfully repro- 
pape As in — & Empe 1; but EmpE eliminated the errors, pointed out by TALLovist, 
in editions 2- 


34. GREENHILL’s first zero of J_1/3(x). 


In Cambridge Phil. So., Proc., v. 4, 1881, p. 68, G. GREENHILL gives the first positive 
root of J_3(x) = 0 as 1.88; this was reproduced in the first edition (1909) of JAHNKE and 
Empe, Tables of Functions, p. 106, and in A. Dryntx, (a) Russkoe fiziko-khimicheskoe 
Obshchestvo, Zhurnal, Chast’ fizicheskaia v. 42, 1911, p. 4 of “Tablitsy funktsii Besselia 
Ja i J44”; (b) Archiv Math. Phys., s. 3, v. 18, 1911, p. 338. 

J. R. Atrey gave the more accurate value 1.8663, in Phil. Mag., s. 6, v. 41, 1921, p. 203. 
Y. IkepA in Z. angew. Math. u. Mech., v. 5, 1925, p. 81 derived the value 1.86635 0858. 
In a table of the first five zeros of Ja} and Ji3, Dinnik gave the value 1.860 (Akad. Nauk, 
U.R.S.R., Kiev, Prirodnicho-tekhnichnii viddil, 1932, p. 12). D. H. L. writes, ‘‘more ac- 
curately this zero is 1.86635 08588 73895 . 

R. C. A. 


35. K. Hayasut, A. Sieben- und mehrstellige Tafeln der Kreis- und Hyperbel- 
funktionen und deren Produkte sowie der Gammafunktion . . . , Berlin, 
1926, table of sin-' x for radian argument x = [0.00000(0.00001)0.001; 
20D], [0.0010(0.0001)0.0999; 10D] and [0.100(0.001)1; 7D], p. 5-86; 
B. Fiinfstellige Funktionentafeln . . . , Berlin, 1930, table of sin x, for 
x = [0.00(0.01)1; 5D], p. 2-4. 


With the exception of the 100 values given to 20D all of Hayashi’s inverse sines were 
proof-read against the Project’s manuscript of the same function, see UMT 7, p. 94-95, 
and the following errors were discovered: 
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A 
x for read = for read 

0132 .01320 03883 01320 03834 479 4995154 4995152 
.0168 .01690 07904 .01680 07904 A8i 5017937 5017950 
.0257 .02570 28099 .02570 28299 5562452 5562438 
0532 05322 51278 05322 51268 543 5744904 5740056 
.0558 .05585 89975 .05582 89975 548 5798714 5799714 
0582 05823 29154 .05823 29064 .6298207 .6298209 
.0626 .06264 09575 .06264 09580 .640 .6944985 6944983 
0627 .06276 11548 .06274 11548 .657 7171234 7168325 
.0660 .06604 80057 .06604 80102 .676 -7422209 7423209 
0691 06915 51094 06915 51084 683 -7519620 -7518620 
0693 .06935 56180 06935 55890 .690 .7614939 7614891 
0732 .07326 53661 .07326 55287 .696 -7698129 -7698116 
.0759 .07597 29638 .07597 30638 -705 .7824131 .7824231 
.08871 63294 08871 63291 .707 .7852467 7852472 
.0891 .08921 83127 .08921 83145 .732 8212534 8212529 
.132 .1323913 .1323864 744 8390418 8390370 
1525922 1525915 .776 .8872990 8882990 
.198 1993227 1993171 837 9922139 .9917776 
.267 .2702835 .2702787 856 1.0274823 1.0274821 
.387 3973754 3973758 1.0411481 1.0411781 
.390 3906316 4006316 .914 1.1530364 1.1530362 
4202672 4202624 927 1.1863332 1.1863334 
457 4746204 4746195 .933 1.2026605 1.2026609 
467 4858971 4858950 .948 1.2473617 1.2468920 
470 4891804 4892908 951 1.2564570 1.2564542 
472 4915880 4915580 


In addition to the above, errors of a unit in the last place (a) a unit too large, (b) a unit 
too small (including doubtful rounding off), occur at the following values of x: 

(a) x=.0015, .0109, .0134, .0138, .0177, .0238, .0381, .0389, .0402, .0539, .0647, .0666, 
.0746-.0748, .0817, .0819, .0924, .0963-.0969, .304, .329, .359, .422, .468, .506, .531, .539, 
-567, .595, .628, .691, .694, .703 and .742. 

(b) x=.0209, .0272-.0279, .0305, .0512, .0614, .0616, .0619, .0707, .0715, .0836, .0876, 
.0878, .0889, .0903, .0914, .0916, .0923, .0927-.0929, .0947, .0949, .0957, .0973, .0976, .0980, 
.0981, .0985, .0989, .167, .180, .181, .195, .346, .368, .369, .395, .453, .458, .516, .520, .529, 
.579, .582, .583, .602, .617, .633, .635, .637, .667, .740, .773, .986 and .994. 

B 


x=.21, for .21156, read .21157; x=.39, for .39063, read .40063; x=.47, for .48918, 
read .48929. 


Although in 1932 Hayashi issued a pamphlet of corrections to his Sieben- und mehr- 
stellige Tafeln he failed to note any errors in his sin x values. 


H. E. Sauzer and E. Isaacson 
Mathematical Tables Project 
New York City 


EprtoriaL Notes.—It seems desirable to elaborate the incidental reference to the 12-page 
pamphlet in which Hayashi published an extensive list of corrections. It is entitled Be- 
richtigung in Hayashis sieben |sic] u. mehrstellige Tafeln (1926), Fukuoka, 1932. 19.1 x 26.4 

he photoprint of this pamphlet in the Library of Brown University was made from 
an original in the Library of Columbia University. There are 8 pages (with 3 columns on 
a page) of corrections; in e* and e~*, 1} cols. (about 75 corrections); in sin x and cos x, 
about 11 cols. (550 corrections), etc. Eleven corrections of the Berichtigung, and six addi- 
tional corrections have been made by hand (presumably Hayashi’s). L.J. c. writes as follows, 
however: ‘‘ This pamphlet has the same percentage of error as the original volume and, what 
is more, there are some cases where things that were right in the original are put wrong by 
these so-called corrections.” 

We may also add notes cupatns errors in radian tables of tan“ x in the above men- 
tioned works of Hayashi, A and B. In A the tables are for x = [0.00000(0.00001)0.001; 20D], 
[0.0000(0.0001)0.0999; 10D], [0.000(0.001)2.999; 8D], [3.00(0.01) 10.0(0.1)20(1)50; 7D]. In 
his Tables of tan x and log (1+-x*), (Tracts for a no. XXIII, 1938), L. J. Comrie 
states that Hayashi’s tables were first checked through the first quadrant, i.e. up to 
x= 1.570, and 100 errors exclusive of errors in the last decimal were found. Comrie remarked 
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further that over 100 errors were found for later values of x. In the above mentioned Hayashi 
pamphlet 71 corrections are given for these tables. In B, 13 errors not mentioned in his 
errata occur in the table for tan! x, x = [0(0.01) 10.00; 5D], and are listed in The Table of 
Arc Tan x by the MATHEMATICAL TABLES PROJECT (New York, 1942); see RMT 90. 


36. P. R. E. JAHNKE and F. Empe, Table of Functions, first-fifth eds., 
1909-1943. Compare RMT 113; MTE 21, 23, 33, 34, 37. Page refer- 
ences are to the first and fifth editions. 


P. 54 (62), F(15°, 32°), for 0.5604, read 0.5603. P. 56, F(60°, 5°), for 0.08745, read 
0.08735, error in first ed. only. P. 61 (69), (15°, 16°), for 0.2788, read 0.2790; E(20°, 16°), 
for 0.2786, read 0.2788; E(25°, 16°), for 0.2790, read 0.2786. P. 64 (72), E(50°, 64°), for 
1.0072, read 1.0007; E(80°, 64°), for 0.9027, read 0.9072; E(90°, 86°), for 0.9926, read 0.9976. 
P. 68 (85), E(8°), for 1.5630, read 1.5632, error in first and second eds. only; K(86°48’), 
for 4.2744, read 4.2746; K(87°36’), for 4.5619, read 4.5609; K(89°36’) for 6.3504, read 6.3509. 
Except in the two cases noted these errors occur in all five editions. 

S. P. GLazENaP, Matematicheskie i Astronomicheskie Tabliisy, Leningrad, 1932, p. 214- 
215. 


37. P. R. E. JAHNKE and F. Empe, Table of Functions, third, fourth and 
fifth eds., 1938-1943. Compare RMT 113, and MTE 36. 


On p. 79 the following recurrence formula is given for the complete Legendre elliptic 
integral: 

4(n + 1)? fEx"dx — (2n + 3)(2n + 5) = 2x"4[E — (2n + 3)(i — x)K]. 
It is readily verifiable that the right hand member should be 


2x"*1{[2n + 1 — (2m + 3)xJE + (1 — x)K}. 
I. OPATOWSKI 
Illinois Institute of Technology 


Epitor1aL Note.—This formula was given correctly in the source indicated by Empg, 
namely: K. F. MUb.uer, Archiv f. Elekirotechnik, v. 17, 1926, p. 337. 


38. MATHEMATICAL TABLES PROJECT, Tables of Arc Tan x, New York, 1942. 
Compare RMT 90. 


P. 27, x = 2.554, 2.557-2.559, insert missing decimal points in function. 
P. 27, x = 2.570, insert missing decimal point in argument. 

P. 41, x = 3.999, insert missing decimal point in function. 

P. 68, x = 6.645, insert missing 1 in function. 

P. 160 (some volumes, not all) x = 5450, insert missing 1 in function. 


Mrs. RutH CAPUANO 
50 Church St., New York City 


39. MATHEMATICAL TABLES Project, Tables of the Exponential Function e?, 
New York, 1939. 


The integral part of the entry corresponding to the argument 2.3979 on p. 240 should 
be 11 in lieu of 10. 
A. N. Lowan 


40. MATHEMATICAL TABLES Project, Table of Natural Logarithms, 1941. 


V. 1, p. 82, argument, for 9023, read 8023; v. 2, p. 355, argument, for 85201, read 85301. 
J. W. Jr. 
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41. A. J. Tuompson, Tables of the Coefficients of Everett's Central-Difference 
Interpolation Formula, London, second ed., 1943. See RMT 148. 


An error in the first edition (1921) was pointed out by CHAPPELL (see MTE 31). 

Page 10, e2 for ¢ = 0.468, for .06091 51280, read .06091 61280. A recent comparison 
by Mr. E. S. Davis has shown that there are no other discrepancies between the two 
editions. The columns ¢: and « have been checked by Miss D. P. K1LNer by comparison 
with the appropriate columns of my Burroughs-script tables of 4-point and 6-point Lagran- 
gean coefficients. From these checks, and Thompson’s known high standard (no errors are 
known in his monumental Logarithmetica Britannica), there is every reason to believe that 
the rest of his figures are correct. 

The following trivial corrections in the text have been pointed out to me in letters from 
Thompson and J. C. P. MIL_er. 


Page line for read authority 
Viii 12 J.C.P.M. 
32 —2 to € to 


42. H. WEBER, Theorie der Abelschen Functionen vom Geschlecht 3, Berlin, 
1876. The following corrections in Table II, p. 183 (The complete sys- 
tem of odd characteristics), were given in a letter, dated 10 February 
1933, from the late H. S. WuiTE to the late W. F. OsGoon. 


“Left-hand column, 7th characteristic from the top, instead of 
should read (at least this is one of the 8 correct systems) 
Bi Bz Ba Bs Bs Bs Br 


my 111 001 110 101 010 100 = «O11 
111 111 011 100 011 111 100 110 


Otherwise, replace this characteristic, no. 7, left, by the no. 11, right, inverted, with its 
76’s also inverted. 
100 


“In system 9, right-hand column, replace & by 110° 


In system 8, right-hand column, replace by 


In system 12, right-hand column, replace by 


In system 13, right-hand column, replace fs by al 
010 
In system 11, left-hand column, replace 8, by 1" 


In system 17, right-hand column, replace & by ~d 


It is singular that Weber should have let so many errors slip by! And notice that five of 
them are in column #,. That must indicate something about his field of vision.” 


43. H. E. H. Wrincu and D. M. Wrinca, Phil. Mag., s. 6, v. 47, 1924, 
p. 63, Io(30) = 7.81674 X 10"; R. C. CoLWELL and H. C. Harpy, Phil. 
Mag., s. 7, v. 24, 1937, p. 1046, Io(30) = 5.70 X 10". 


The Epirors brought these contradictory results in tables, already referred to in 
MTAC (p. 139-140), to my attention, and suggested that it would be a matter of interest 
to have the correct result determined. 
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In the asymptotic expansion 


upon setting x = 30, and evaluating 25 terms of the series, each to 22 decimal places, I 
found the sum of the terms in the braces to be 


1.00424 76530 20713 59155(8). 


This calculation was performed twice,—first by means of a calculating machine, the 
second time long-hand. The results agreed perfectly. 

The value of e” was found in the W.P.A. Tables of the Exponential Function e (1939), 
p. 533, to 19S. Subsequently I checked and extended this approximation to 38S. I have 
unpublished values of x# and 1/24 to 317 and 310 decimal places respectively. These were 
calculated from x and 1/z, respectively, and the product of the square roots was formed 
with the assistance of a machine, and was found to consist of a sequence of 309 con- 
secutive 9’s. Consequently I have great confidence in the accuracy of these roots. 
The value of 154 was determined to 40D. By multiplication I found ¢°154/3074 
= 7.78366 06884 04464 04193 55906 75 X 10"; and finally, (30) = 7.81672 29782 39774 8972 
X 10", correct to 20S. 

Hence the value of Wrinch and Wrinch is nearly correct, the error being in the last 
significant figure which they give. The value of Colwell and Hardy is entirely incorrect. 

J. W. WreEnca, Jr. 
4211 Second St., N.W. 
Washington 16, D. C. 


UNPUBLISHED MATHEMATICAL TABLES 


References have been made to Unpublished Mathematical Tables in 
RMT 138 (Davis, MILts), and RMT 157 (Davis, Matu. TABLEs Project). 


21[B].—Tables of Fractional Powers, Mss. prepared by, and in possession of, 
the MATHEMATICAL TABLES PRojEcT, 50 Church St., New York City. 


Of these 12 tables 6 are of A*, and 6 of x*, as follows: 


I: A = 2(1)9, and x = [0.001(0.001)0.01(0.01)0.99; 15D]. 
II-IV: A = 10, z, and Euler constant, x = [0.001(.001)1.000; 15D]. 
V: A = 0.01(0.01)0.99, and x = [0.001(0.001)0.01(0.01)0.99; 15D}. 
VI: A = p10, » = the primes between 101 and 997, and x = [0.001(0.001)0.01- 
(0.01)0.99; 15D]. 
VII-XI: x = [0.01(0.01)9.99; 15D] and a = +}, + }, + 3, + 3, + 3. 

XII: a = 0.01(0.01)0.99. For any a the values of x* were computed at the interval 
0.001 in x approximately up to that value of x for which the derivative of the 
function (which is © at x = 0) is in the neighbourhood of unity. All entries are 
to 7D. 

A. N. Lowan 


22[D, H].—Roots of the equation tanx + ax = 0, a = 3/(8r), 3/(12m), 
3/(167). Ms. in possession of the Department of Aeronautical Engi- 
neering, University of Michigan, Ann Arbor, Mich. Compare Q. 8. 


The first three roots of each of these equations were calculated to 5S, and the next 
five roots to 3S. First, we obtained an approximate value x,°, graphically, and then applied 
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a correction 
én =— (tan x,° + ax,°)/(sec* x,° + a). 
The corrected values x,° + ¢, are about as accurate as we need them. 


Orto Laporte 
University of Michigan 


23[L].—Hansrajy Gupta, Table of Liouville’s function and its sum. Ms. in 
possession of the author, Government College, Hoshiarpur, Punjab, 
India. 


Liouville’s function \(m) may be defined for the positive integer argument m = p:% 
2° --+ p** to be + 1 or — 1 according as m = a1 + a2 +---+ a: is even or odd. This 
table gives \(m) and also the sum 


L(n) = M1) +---+ A(n) 


for n < 20 000. The function L(n) is clearly the excess of the number of those integers n 
which have an even number of prime factors over the number of integers having an odd 
number of prime factors. The reason for computing this table is to verify or refute an im- 
portant conjecture of Pétya (Deutsche Math.-Ver. Jahresb., v. 28, 1919, p. 38) to the 
effect that this ‘‘excess” is really negative or zero for n > 1. The table shows that the 
conjecture is correct as far as nm = 20 000 but does not show that the conjecture has safely 
passed its worst trials. In fact one finds that L(48 512) = — 2. Data, taken from this table, 
on the behavior of L(m) have been published by the author, Indian Acad. Sci., Proc., v. 
12A, 1940, p. 407-409. Current interest in Pélya’s conjecture (which implies the Riemann 
Hypothesis) has been heightened by the recent results of A. E. INGHAM, Amer. J. Math., 
v. 64, 1942, p. 313-319. The corresponding problem for the companion function y(n) of 
MO6stvus and its sum is also unsolved and has been the subject of much tabulation; see 
Guide to Tables in the Theory of Numbers, Nat. Res. Council, Bull. no. 105, Washington, 
D. C., 1941, p. 9. 
D. H. L. 


NOTES 


17. DESIRABLE MATHEMATICAL TABLES FOR REPUBLICATION.—The 
Office of Alien Property Custodian has licensed, during the past several 
months, the reprinting of scientific and technical books, of enemy origin, 
which are not available in a quantity sufficient to meet the demands of the 
wartime operations of science and industry. In RMT 79, 113, 128, references 
have been made to three volumes of this kind. Before definite decision can 
be made regarding the licensing of additional Mathematical Tables for 
republication, it is necessary for the Custodian to be informed about the 
extent of the need of such tables and to receive suggestions of specific titles 
for consideration. This can be accomplished if suggestions of specific sig- 
nificant tables, or any inquiries, are sent by individuals to the undersigned 
at Division of Patent Administration, Office of Alien Property Custodian, 
Washington, D. C. 

H. H. SarGEAnt, Chief 


18. Phil. Mag. TaBLEs, SuppL. 1 (see MTAC, p. 135-141).—In H. C. 
PLuMMER, “The numerical solution of a type of equation,” s. 7, v. 32, 1941, 
p. 505-512, roots are found of the following transcendental equations, 
among others, of the type tan x = xf(x), where f(x) is a single-valued 
function: 
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I. For tan x = x/(1 — x*), which relates to electrical waves on a sphere 
(J. J. THomson, Notes on Recent Researches in Electricity and Magnetism, 
Oxford, 1893, p. 373), are given the first 5 positive roots to 4D. The first 3 
of these roots agree with those of L. SILBERSTEIN, Bell’s Mathematical 
Tables, London, Bell, 1922, p. 97. Silberstein’s fourth root is appreciably 
in error. Before either of these results had been published, H. Lams gave 
the first 6 positive roots to 4D in London Math. So., Proc., v. 13, p. 202, 
1882. 

II. tan x = 2x/(2 — x*), an equation found in connection with discus- 
sion of sound waves, or the free oscillations of gas, in a rigid spherical shell 
(RAYLEIGH, The Theory of Sound, v. 2, 1878, p. 232; second ed., 1896, p. 
265). Five positive roots to 4D are given. For this equation also LAMB gave 
(l.c.) 6 positive roots to 4D. 

III. Another equation considered by Rayleigh (tbid., v. 2, 1896, p. 266) 
is tan x = x(x? — 9)/(4x* — 9) for which he gives one non-zero solution to 
4D; Plummer finds four positive roots to 4D. 

IV. The equations cos x-cosh x = + 1, occur in connection with the 
lateral vibration of bars (RAYLEIGH, ibid., v. 1, 1877, p. 222-224; second 
ed., 1896, p. 278). Four positive solutions, to 6D, are found for each of these 
equations. In this connection Rayleigh’s errors in the first edition are cor- 
rected in the second. The oversight was noted by A. G. GREENHILL (Mess. 
Math., v. 16, 1886, p. 119) who has given another interesting application 
of the same equations. More roots and general formulae are given in E. P. 
ApAms, Smithsonian Mathematical Formulae, Washington, D. C., 1939 


(1922), p. 86. H. B. and R.C. A. 


19. Phil. Mag. TaBLes, SuppL. 2.—In H. Carrinoton, “The fre- 
quencies of vibration of flat circular plates fixed at the circumference,” 
s. 6, v. 50, 1925, p. 1261-1262, the 16 roots, < 16, of Jn4i(x)-In(x) 
+ In41(x)-Jn(x) = 0, for n = 0, 1, 2, 3, are given to 5S. The first ten roots 
for each of these values of m are given by J. R. Arrey, ‘‘The vibrations of 
circular plates and their relation to Bessel functions,” Phys. So. London, 
Proc., v. 23, 1911, p. 227; for m = 0 these roots are given to 4D, but for 
the others, to 3D. Carrington’s first root for » = 0 was 3.1961; according 
to Airey this should be 3.1955. Of the 10 roots given in JAHNKE and EMDE 
(1933), p. 283, and 1938-43 editions, p. 234, 8 are erroneous to the extent 
of 1 to 20 points in the last decimal places. RCA 


20. A Root OF THE EQUATION 10 log x = x.—In Mathematical Gazette, 
v. 15, p. 367, 1931, C. V. Boys gives what purports to be a 60D approxima- 
tion to the irrational root of this equation. I checked this value and found 
it correct to only 38D. By Newton’s method I deduced a value of x to 66D 
and then in a second calculation determined the logarithm of this last 
approximation. It developed that my result was correct to 65D. Inci- 
dentally I employed the rational approximation x = 27.5.7/3*.112. The 
correct value for x, to 65D, is as follows: 


1.37128 85742 38623 53686 13621 06299 68995 88428 54404 84225- 
70704 08723 85385. 
J. W. WRENCH, Jr. 
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Epitortat Note.—Part of this value is the first in a table of ten “coincidental log- 
arithms”’ listed by D. M. KNappPeEN,' solutions (8 of them irrational) of the equations 
10" log x =x, n=1, 2, ---, 10. To 50D it agrees with the result of Mr. Wrench; in the re- 
maining three places there is disagreement, 650 instead of 707. For n=8 the value given 
for x, to 66D, is as follows: 


8.95191 59982 67846 23159 99873 77688 49072- 
72932 33707 59495 47848 84326 51391 4 Xx 108. 


(Upon bringing this result to the attention of Mr. Wrench he verified its accuracy through 
64D and then added the following 16 digits: 0 73919 87374 90256.] Knappen states that 
these remarkable numbers were called “constants’’ by OLIVER ByRNE, = certainly dis- 
cussed such numbers in at least two publications,” of 1849 and 1864. In each of these places 
Byrne considered the problem “‘to find a number whose common logarithm is composed of 
the same digits and in the same order as itself.’’ In the first he gives the ten values, each to 
15D, but the last two or three figures of eight of these are incorrect. Similar remarks apply 
to the seven values given in the second publication, where Byrne recklessly writes “‘no 
known development but the dual will establish” them. As yet we have not been able to 
connect with Byrne, the 53 to 66D values given by Knappen. Perhaps some reader can 
supply this information. 


1D. M. Knapren, Math. Mag., Washington, D. C., v. 1, p. 202—the number con- 
aes | this page is dated Oct. 1884, although published in August 1887. 

20. Byrne, (a) Practical, Short, and Direct Method of Calculating the Logarithm of any 
iven Number, and the Number corresponding to any given Logarithm. London and New 
ork, 1849, p. viii, 14-15, 52, 65-76. Compare J. encase, Bibliotheca Tabularum 

Mathematicarum, Sy 1926, p. 121-122. 
(b) Dual Arithmetic. A New Art... . New issue with a Complete Analysis, London, 
1864, Analysis, p. 75-78. 


QUERIES 


8. RoOTs OF THE EQUATION tan x = cx.—As early as 1748 in v. 2 of his 
Introductio in Analysin Infinitorum, Lausanne, p. 318-320, EULER solved 
the problem “‘Invenire omnes Arcus, qui Tangentibus suis sint aequales,”’ 
by deriving a general formula, from which he found the first 10 non-negative 
roots of tan x = x(c = 1). The second of these roots was given in the form 
“3.90°— 12° 32’ 48” ” [4.49340834 instead of 4.49340964], and the tenth as 
“*19.90°— 1° 55’ 16”. These results may be regarded as correct to 5D. 
Js2(x) = 0 is equivalent to Euler’s equation tan x = x. In 1886 LoMMEL 
listed the first 17 non-negative roots, to 6D (Akad. d. Wiss., Munich, 
Math.-naturw. Abt., Abhandlungen, v. 15, p. 651). This list, abbreviated to 
4D, was published in the first, second and fifth editions of JAHNKE and 
EmDE (see RMT 113), and in E. P. ApAms, Smithsonian Mathematical For- 
mulae . . . , Washington, D. C., 1939 (1922), p. 84. The first 36 roots, to 
4D, are given in K. Hayasni, Fiinfstellige Funktionentafeln, Berlin, 1930, 
p. 52. The 18th to the 36th roots, to 6D, are given in L. SILBERSTEIN, Bell’s 
Mathematical Tables, London, Bell, 1922 (= Synopsis of Applicable Mathe- 
matics with Tables, New York, Van Nostrand, 1923), p. 97. In 1827 Caucny 
found the first 4 roots to 7D when c = 1, the first 5, to 6D, when c = 4, 
and the first 4, to 4-6D, when c = 8/5 (‘Théorie de la propagation des ondes 
a la surface d’un fluide pesant d’une profondeur indéfinie,” Acad. d. Sci., 
Paris, Mémoires prés. par divers Savants, s. 2, v. 1, p. 198-209; Oeuvres, 
s. 1, v. 1, Paris, 1882, p. 204-212). Compare UMT 22. What other solutions 
of tan x = cx are known? 

H. B. and R. C. A. 


9. System oF LINEAR EQuaTions.—Where may one find an adequate 
treatment of the method of Gauss and SEIDEL for solving a system of n 
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linear equations in m unknowns by successive approximations? The dis- 
cussion given in WHITTAKER and Rosinson, The Calculus of Observations 
(London, 1924, and third ed., 1940, p. 255-256), is not satisfactory. The 
part purporting to show that the process always improves a trial solution 
suffers the following simple exception: 


2a+y=1, x+3y=-1. 


Here the initial solution x = 1/2, y = — 1/3 is not improved by replacing 
x by 2/3 as required by the process. 
D. H. L. 


QUERIES—REPLIES 


8. TABLEs or N*?(Q5, p. 131).—Another table for three-halves powers 
of numbers to more than three places is T. 70, p. 290 of J. T. FANNING, A 
Practical Treatise on Hydraulic and Water-Supply Emergency, tenth ed., 
New York, 1892, where N = [0.04(0.01)0.20(0.02)1.0(0.1)4; 4D]. 
H. B. 


CORRIGENDA 


P. 2, 1. 31, for Reply to Query 6, read Reply-to-Query 6. P. 6, 1. 6, for v. 4, read v. 14. 
P. 9, 76 for CHAPMAN, read Cuapin. P. 14, |. 5 from bottom, for 0.001, read 0.0001. 
P. 15, 1. 6, add Also, p. 224-224c, sin x, cos x to 10D, log sin x, log cos x to 5D, 
x = 0(.1)10, 0(1)100. P. 15-16, omit references to HAyAsui tables of sin 4x2, cos }xz, 
1. 13-14 from bottom of p. 15; also to tables of KOLKMEIJER and BUERGER, top of p. 16. 
P. 16, 1. 8 from bottom, for Spoon, read Spon. P. 18. 1. 1 and 2 from bottom, for 6D, 
read 6D-7D. P. 19, 1. 3 from bottom, for x, ---3D], read x = [0.00(0.01)1.0(0.1)10(1)- 
100(10)1000; 3D]. P. 20, footnote, 1. 6, after ‘‘109.’’ insert With the aid of the entries 
presented the logarithms of all numbers N = 1(1)109 are readily found. P. 47, 90, 1. 3, 
for State, read City. P. 69, 2, 1. 3, for Houghton, read Haughton; 3, |. 1, for 12S, read 
10S-12S; 5 and 6, for with differences, read with first differences. 

P. 70, 8, 1. 2, for 10D, read 9D-10D; 1. 4, for 0(1/2)(13/2), read 0(4)63; 1. 5, for in read 
3nzx, [this was a mistake in the Report]; 10, 1. 3, for by degrees, read at three-degree 
intervals; 12, |. 3, for 80°1, read 80°.1. P. 73, 44, 1. 2, Ei in roman, not ital.; 49, 1. 4, 
for 0.0(0.1)10.0 read 0.0(0.1)7(1)10. P. 74, 52, 1. 20, for and read and 
56, 1. 4, for 120, read 12.0. P. 96, in UMT 9, totals, make the following changes: 390 
for 391; Poulet 65 (for 68); Escott 233 (for 235); and add Poulet and Gérardin 4 (1929). 
P. 109, 1. 17-18, for Ji(17), read Ji(x:7); 1. 20-22, for these lines read, the roots of 
Ji(x) Ni(kx) — Ji(kx)Ni(x) = 0 on p. 204 of nos. 3-5, p. 274 of no. 2, and p. 162 of no. 
1, the first three roots for the value k = 2 should be 3.1917, 6.3116, and 9.4446 accord- 
ing to values given in Muskar, ---. P. 108, 1. 17, for Debeye, read Debye. 

P. 125, 1. 20-23, for numbers, read figures. P. 138, 26, 1. 4, for J44(x), read J+4(x); for uncer- 
tain fourth, read approximate fifth; 1. 5, for J43(x), read J+3(x); 1. 5-6, for uncertainties, 
read approximate fifths; 1. 7, for }(m + 1), read }/(m + 1). P. 140, no. 38, for dx, and 
dx?, read dv and dr”. P. 143, 1. 4 from bottom, for einen, read einem. P. 145, for line 8, 
read: place tables for A with D = 0.0000(0.0001)2.000(0.001)4.00(0.01)6.94; and for S 
with D = 0.3000(0.0001)2.000(0.001)4.00(0.01)6.94. P. 157, 1. 16-17, for B,“™(0) and 
B,™ (1), read B,™(0)/n! and B,™(1)/n! P. 161, 1. 11, delete ““P. 54, F(35°, 30°), for 
0.6220, read 6200.’’ P. 161, 1. 13, for 1035, read 1037. P. 164, 1. 11 from bottom, eliminate 
the second ‘‘10;’’. P. 168, 1. 26, for Kiister, read Kiistner. P. 169, 1. 27, read Physical; 
1. 6 from bottom, for kkadratov, read kvadraiov; |. 4 from bottom, read Izdatyel’stvo. 
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EDITORIAL NOTICES 


The addresses of all contributors to each issue of MTAC are given in 
that issue, those of the Committee being on cover 2. The use of initials only 
indicates a member of the executive committee. 


The next number of MTAC will be a greatly enlarged special issue 
wholly devoted to a Guide to Tables in Bessel Functions. 


Volume I is to contain at least 8 numbers continuously paged. 
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